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A general method to derive the diagonal representation for a generic matrix valued 
I quantum Hamiltonian is proposed. In this approach new mathematical objects like 

non-commuting operators evolving with the Planck constant promoted as a running 
Q^i variable are introduced. This method leads to a formal compact expression for the 
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INTRODUCTION 



diagonal Hamiltonian which can be expanded in a power series of the Planck con- 
stant. In particular, we provide an explicit expression for the diagonal representation 
of a generic Hamiltonian to the second order in the Planck constant. This last result 

in ; 

\ is applied, as a physical illustration, to Dirac electrons and neutrinos in external 
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■ The spacetime evolution of some important quantum systems is governed by a multicom- 

ponent Schrodinger like equation whose Hamiltonian is a matrix valued operators. In many 
cases, these systems are too complicated to be solved directly using the Schrodinger equation 
and some simplifications are required. For systems displaying a separation of scales in terms 
of slow and fast degrees of freedom a major simplification is to treat quantally the fast ones 
while the slow ones can be approximated classically or semiclassically. For a slow cyclic 
variation of the classical part. Berry has shown that the wave function of the quantum part 
gets an additional geometric phase factor [1]. It has also been understood, in the context 
of the Born-Oppenheimer theory of molecules [2] [3] [4] that the geometric phases induced 
by the quantum part of the system provide a correction, reacting on the classical part with 
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geometric Lorentz and electric forces [5] [6]. Going beyond this first correction has proved 
extraordinarily difficult due to the intricate entanglement of noncommuting operators [7] . 

In fact, many physical systems share a very similar mathematical structure with molec- 
ular systems. For instance, the translational movement of a Dirac electron in a slowly 
varying external field can be separated from the spin evolution and therefore considered 
semiclassically. 

In this paper, we adress the question of the dynamics governing the slow part of an 
arbitrary quantum system and present a formalism for calculating the Hamiltonian of this 
slow part to all orders in the Planck's constant. Of course, we do not expect the conver- 
gence of the scries in general, since for most quantum systems the limit h is singular, 
which means that an exact quantum behavior can not be approximated even by an infinite 
expansion. 

The philosophy behind the approach developed here, consists in a mapping of the initial 
quantum system to a classical one which can be easily diagonalized and then to return to the 
full quantum system. The method requires the introduction of new mathematical objects 
like non-commuting operators evolving with the Planck constant promoted as a running 
variable. This leads to a differential calculus on a non- commutative space which shows some 
similarities with the stochastic calculus. This innovation is the clue to derive the formal 
exact diagonal representation for any arbitrary matrix valued Hamiltonian. 

The results presented here extend the work presented in [8] , where a differential equation 
with respect to h was proposed for the required diagonal Hamiltonian. The semiclassical 
limit, which is often enough to get physical insight for the problem considered, was obtained 
by a straightforward integration of this differential equation and a general semiclassical 
formula for an arbitrary diagonal Hamiltonian in terms of covariant operators and commu- 
tators between Berry connections was obtained. This formalism allowed to consider several 
straightforward applications. These ones included the deffection of light by the interaction 
between its polarization state and external inhomogeneities, such as a varying refractive 
index or a gravitational field [9]; Dirac electrons and Bloch electrons in crystals [10], leading 
to new effective geometrical forces and a clarification of the Peierls substitution [11] as well 
as to a generalization of the Bohr-Sommerfeld quantification rule [12]. In principle, one 
could solve the differential equation to the desired order in a series of h, but although this 
was done for simple examples to order (see [8]), we could not derive such an expression 
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for an arbitrary Hamiltonian. 

The main advantage of the new method presented in this paper is that being not based 
on a differential equation, it directly leads to a formal exact diagonal representation for any 
arbitrary matrix valued Hamiltonian, which, as proven in this paper, is an exact solution of 
the differential equation of [8]. It appears also, as one could expect, that the expansion in 
series of H for the exact diagonal Hamiltonian derived here is much more easier to obtain 
than by the successive integration of the differential equation of [8] . In particular, we provide 
an explicit formula for the diagonal Hamiltonian to order that can be directly applied in 
principle to any physical systems. 

In this paper, for the sake of clarity and to keep the length reasonable, we will consider 
only two simple but still physically relevant examples. More applications can be found in 
[13], where interacting Bloch electrons are considered, and in [14] for a generalization of the 
Born-Oppenheimer approximation. 

It is worth mentioning that some other diagonalization procedures exist, each of them 
having its range of validity, advantages and defects. The paradigmic example is provided by 
the Foldy-Wouthuysen (FW) representation of the Dirac Hamiltonian for relativistic parti- 
cles interacting with an external electromagnetic field. In this representation the positive and 
negative energy states are separately represented and the non-relativistic Pauli-Hamiltonian 
is obtained [15]. Actually even if several exact FW transformations have been found for 
some definite classes of potentials [16] [17] [18], the diagonahzation is a difficult mathematical 
problem requiring approximations, essentially a perturbation expansion in weak fields. A 
generahsation of the FW transformation developed by Blount and also based on an expan- 
sion in weak electromagnetic fields, is in principle applicable to Bloch electrons and Dirac 
electrons [19]. But this method is limited by construction to weak external perturbations 
and provides a formal series expansion which leads to cumbersome calculations for practical 
applications. Recently a variant of the Foldy Wouthuysen transformation valid for strong 
fields and based also on an expansion in h of the Dirac Hamiltonian was presented [20] . 

Among other approaches, Weigert and Littlejohn developed a systematic method to diag- 
onalize general quantum Hamiltonian in a series expansion in h [7]. It leads also to a formal 
series expansion written in terms of symbols of operators which also makes the method 
complicated for practical applications. In particular this method does not provide explicit 
formulas, but instead is an algorithm to compute case by case the series order by order. 
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To our knowledge, this approach was only used for Born-Oppenheimer molecular systems. 
Note also, that a different method based on a adiabatic perturbation theory was proposed 
in [21] (and references therein) and applied to the Dirac electron and electrons in a cristal 

at the first order. 

Our approach, although similar in spirit to these methods, is technically very different 
and allows us for the first time to derive the expression of the diagonal representation at 
order of an arbitrary quantum Hamiltonian. This efficiency could be explained by the fact 
that, unhke other approaches, no Weyl calculus and Moyal Algebra are required. Another 
interesting feature of our approach is that it confirms the fundamental role played by Berry 
curvatures since it leads to an effective diagonal Hamiltonian with Berry phase corrections 
as well as noncommutative (Berry connections dependent) coordinates and momentum co- 
variant operators as in previous approaches [8] [10] (see also [22] for Dirac electrons in an 
electric field and [23] for the extension to the full electromagnetic field). The resulting 
generic equations of motion are also corrected by Berry curvatures terms. 

The paper is organized as follows. In the next section we develop the differential calculus 
in noncommutative space. We then derive, in section 3, the diagonalization procedure 
and show the link with the differential equation in section 4. Section 5 is devoted to the 
discussion of the dynamical coordinate and momentum operators. In section 6, we give the 
general diagonal energy operator formula to the second order in h and in section 7, Dirac 
electron in an external electric field and neutrinos in a static symmetric gravitational field 
are considered as an illustration of the general formalism. Last section is for the conclusion. 

II. PRELIMINARIES: 

A. DifFerential calculus on noncommutative space 

Consider a quantum mechanical system whose state space is a tensor product {TZ^)®V 
with V some internal space. In other words, the Hamiltonian can be written as a matrix 
i?(P,R) of size dimy whose elements are operators depending on a couple of canonical 
operators [it!*, P^] — ihSij, the archetype example being the Dirac Hamiltonian with V — C^. 
Our goal is to derive the formal diagonal representation of this matrix valued quantum 
Hamiltonian. By diagonal Hamiltonian it is meant an effective in-band Hamiltonian which 
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has a matrix representation with block-diagonal matrix elements associated with energy 
band subspaces. As in paper [8] the procedure we propose for the removal of the interband 
matrix elements needs the use of some unconventional mathematics compared to the usual 
formalism of quantum mechanics. Since the Planck constant is here considered as a variable, 
operators do not satisfy the usual rules of quantum mechanics. Therefore before directly 
embarking on the diagonalization procedure we first introduce the required mathematical 
tools. For that purpose we start with some definitions and notations. 

1. Running coordinates and momentum operators 

To begin with, we introduce a formal space of non commuting infinitesimal operators 

dX^ = {dRl^,dP^} i = 1,2,3 indexed by a continuous parameter a, that satisfy the 
following infinitesimal Heisenberg algebra (with a reversed sign) : 

[dRl,, dPi] = -ida6a,aAj and [rf^, dlV^] = [dP^, dPi,] = 0. (1) 

Then, we define a set of running coordinate and momentum operators by witting the fol- 
lowing formal sums : 

Ri^R^- TdR,, and P^^ P' - f dP{ (2) 

J a J a 

with the convention di?^ — R^^ — R!'^_^ and dP"^ — Pa — Pa-da- This "downward" choice 
of the differential element notably imphes the commutation rules [i?^, dP^ = [dR\, P^] = 
which turns out to be absolutely necessary to develop later on a differential calculus on this 
noncommutative space. 

For a — hwe recover the usual canonical operators R^ = i?^and = P^ which satisfy 
the canonical Heisenberg algebra [P*, P-'] = ihSij, whereas the running operators satisfy 

[Ra, Pi] = and [Rl, Ri] = [P^, Pi] = 0. (3) 

Note that in this paper we will never consider the algebra of the operators P^ and P^/ for 
a ^ a', which from the definition Eq. (2) is clearly not a Heisenberg one. 

The variables dX^ have to be understood as fictitious variables that make the link between 
quantal operators {a — h) and classical variables (a = 0). As will be explicit later on, their 
role is to transport our quantum system to a formal classical one that can be in general 
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easily diagonalized, and then back from the formal classical one to the required quantized 
system. By writing dX.a — ^Jda/K9ia with a normalized canonical operator we see that 
the infinitesimal quantities rfXo, are actually of order \/'da. However, having this in mind, 
we will never use the notation and always work with dX^. 

2. Differential algebra 

For the sequel, we need to define the differential dF (X^, a) of an arbitrary function 
F(X.a,Oi) where = {i?^,P^}. For this purpose, consider the operators i?^ — dR\, 
and P^^ P'- J^dPl as acting on a space W ^ {V ® L'^ (TZ^)) (g) {®a<hL'^ ('^^)a) which is 
the tensor product of V and an infinite number of copies of (7^^)- The tensor product 
V ®L'^ (7^^) and the space refer respectively to the spaces on which the canonical 

operators {R\P^) and the differential operators dX^^ act. 

Now to be consistent with our convention for the differential element dX^^, the differen- 
tial of the function F(Xq,,q;) in defined in an unusual backward manner as dF(Ka,Oi) = 
F (Xq,, a) — F (X.a-da, Oi — da) . This construction is essential as it permits the commuta- 
tion between the non-differential and the differential elements dX^^. The differential will be 
written as a second order expansion plus some neglected terms: 

dF (X„, a) = VR,FdRi + Vp.FdPl - \ (V«, V«^.F) dRJW, - \ (Vp, Vp^.F) dPldPi 

1 dF 
- 2 (Vii, Vp.F) {dKJPi + dPidR^) + (F (X„ a)) da + —da 

-\- terms of order 3, (4) 

where all expressions in the r.h.s. are evaluated at (Xa,Q;). Note that here, we have kept 
the terms of order square in since they are of order da and thus contribute to the 
differential, whereas higher orders can safely be disregarded as they are negligible when 
the integration over a is considered. These second order terms have been organized in a 
certain form that will be practical for us later. There is nothing to say about the usual 
dRl^dR^^ or dP^dP^ terms. The crossed terms involving products of dP^ and rfi?^ have to 
be taken with care since these two infinitesimal elements do not commute with each other. 
For this purpose we have decomposed the second order terms involving products such as 
dPl^dRP^ and dW^dPl^ in a symmetric part proportional to {dR]^dPi -\- dPidR\) giving the 
contribution — | (Vp^VpjF) [dR^^dPl -\- dPldR]^), and an antisymmetric part proportional 
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to {dR\dP^ — dPldR]^) — —iS^^da corresponding to the bracket {F (X^^, a)) da that we now 
explain in details. The notation (F (X^, a)) (which in [8] was corresponding to the opera- 
tion — lAsi/rnVfl. VpiF (Xa, a)) is defined as a specific procedure on a series expansion of 
F in the variables i?^, in the following way : let F be a sum of monomials of the kind 
Ml (Ra) M2 (Pa) M3 (Rq) .... the Mj being arbitrary monomials in i?„ or P„ alternatively. 
The operator V R^V pj acts on such an expression by deriving all combinations of one mono- 
mial in Rq and one monomial in 'P^- For each of these combinations, insert a dFC^ at the 
place where the derivative Vr. is acting and in a same manner a dP^ at the place where 
the derivative V pj is acting. This leads to an expression with two kind of terms, one kind 
being proportional to the dR\dP^, and the second proportional to dP^dR\. Then rewrite 
this expression in terms of dR^^dP^ + dP^dRl^ and dR^^dP^ — dP^dR^^ = —iS^^da. Then 
{F (Xq,, a)) is defined as minus the sum over i and j of all the terms proportional to da 
in the computation procedure just considered. This definition implies a procedure which is 
clearly dependent of the symmetrization chosen for the expansion of F. 

To make the definition of {F (X.a,a)) clearer, consider some important practical exam- 
ples. If the function F has the following form F = 1 (A (R„) B (P^) + B (P^) A (R„)) 
which corresponds to a frequent choice of symmetrization in R^ and P^, then {F (Xq, a)) — 
J / [Vr^^ (Rq) , VpiB (Pa)]- Another choice of symmetrization leads in general to a dif- 
ferent result. 

For instance, if we rewrite the same function F in a fully symmetrized form in Ra and 
Pa (that is invariant by all permutations in R^ and Pa) which is also often used, we now 
have a different result since (F (X^, a)) = 0. 

Nevertheless, this dependence of {F{X.a,a)) in the symmetrization choice is not aston- 
ishing at all. Actually changing the symmetrization of a function F (Xq, a) introduces 
some explicit terms in a which changes also the term daFda present in the differential Eq. 
(4). As a consequence, neither the partial derivative with respect to a, nor the bracket are 
invariant by a change of form. But, what is invariant is the sum daF + (F). 

We show this assertion by first giving a useful formula for the differential of a product of 
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functions : 

d(F(X,,«)G'(X,,«)) = dFG + FdG - V R,FV R.GdRJR>,-V p,FV p.GdPldPi 

- {Vr^FVp^G) dRJPi - {Vp^FVr^G) dPidK 
= dFG + FdG - V R,FV R^GdRJIV^ - V p^FV p^GdP^dPi 
~ {Vr,FVp.G + VpjFVr,G) {dKJPi + dPidR,) 

-- {Vp^FVr^G - Vp^FVpfi) da (5) 

Now, the invariance of daF + (F) is the consequence of two facts. The first one is 
that the definition of the differential dF(X.a,a) is independent of the way the function 
F (Xq,, a) has been symmetrized. Actually, considering F (X^, a) as a monomial, a change 
of symmetrization amounts to move successively powers of the momentum on the right or 
the left of the position variable. Recursively, one just has to check that moving only one 
power of the momentum does not change the differential. Considering thus a monomial 
of the form F (Xq, a) — MiP^i?^M2 with Mi and M2 arbitrary, move the momentum to 
rewrite F (Xq,, a) — Mi {R^Pl, — ia) M2. Our formula for a product of differential, allows 
to assert that the differentials of the two terms MiP^i?^M2 and Mi {R\Pa — ^tt) M2 differ 
only by the contributions Mirf (P^i?^) M2 and Mid{R\P^ — ia) M2. As a consequence, 
the differential dF (Xq,, a) is independent of the choice of symmetrization if and only if 
d (PiRi) = d {RiP^ - ia). Given that 

d {KP^ - ia) = ajPi + dR^Pi - i (dRJPi + dP^dRi,) + '-da - ida 

and that the dP^,dRl^ commute with P^, P^ the equahty of the two previous lines holds as 
well as the assertion of symmetrization independence. 

This result joined to the second fact that in the series expansion of dF (Xq, a), the vari- 
ables dX^^, dX'^dX^, da, arc independent, and thus that the coefficients of the expansion of 
dF (Xq, a) are uniquely defined, leads directly to the announced result that daF+{F) which 
is the coefficient of da in the expansion of dF (Xq,, a) is independent of any symmetrization 
choice. Therefore it is convenient to introduce a symmetrization invariant derivative Da as 
given by 

P>qP (Xq, a) = daF (Xq, a) + {F (Xq, a)) (6) 
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We end up this paragraph by introducing two formulas that will be convenient later. 
First, we still more compactify our notations for the differential by writing : 

6 ^6 
i=l i,j=i 

+ D^F (X„, a) da (7) 

with i J = 1..6. We also assume that X* = i?^ for i = 1, 2, 3 and = for i = 4, 5, 6. In 
fact, as we will show in the next sections the quantity of real importance for us in Eq. (7) is 
the term proportional to da. Note that in this notation, our previous formula for a product 
of differential Eq. (5) takes a compact form : 

6 

d {F (X„, a) G (X„, a)) = dFG + FdG - ^ Vxj {F) V i^) dX'jXi (8) 

where the last term could of course be developed as before in symmetric and antisymmetric 
parts. 

Second, we also give the bracket formula (.) for a product of two functions. Using the 
procedure defined previously (see also ref. [8]) one obtains the following expression 

< F (X<„ a) G (X„, a) >= {F) G + F {G) - '-Vp.FVr^G + '-Vr^FVp^G. (9) 

This formula shows that the bracket operation for a product F (X^, a) G (X^, a) can also be 
seen as a sort of deformation of the Poisson bracket including some " internal" contributions 
(F(X«,a)) and (^(X^.a)). 

Let us remark ultimately that the term {F (X^^, a)) da in Eq. (4) which is of the second 
order in the derivatives Vij. and Vp< is very reminiscent of the bracket introduced in the 
stochastic calculus. Given the non commutativity of the operators at stake, one should 
in fact rather expect our formalism to be close to the quantum stochastic calculus [25]. 
However, it does not seem at first sight that our objects fit in such a framework which 
deals rather with a formalism of creation, annihilation and conservation operators. A full 
comparison is out of the scope of our paper. 

3. Integration 

Now, we will prove that the formula Eq. (4) allows to express a function F (X;j, h) which 
depends on the physical quantum operators X^ defined at the quantum scale a — h, as the 
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integration of a differential. 
Proposition 

F{-Kn,h)=F{-Xo,0)+ f dF{Xa,a). (10) 

Jo 

Before proving Eq. (10), remark that this expression has the following immediate gener- 
alization : 

F(Xn,/i)=F(X„,a)+ / dF(X^,/3) 

J a 

Proof. 

The proof of formula Eq. (10) is by recursion. Expanding F (Xq,, a) as sum of monomials 
in the canonical variable, it is enough to prove our assertion for a monomial of a certain 
degree in the R\, P^. For monomials of degree one, that is hnear expressions in the R\, 
the result is trivial given our definitions. Now, assume that the result is true for all 
monomial M of, say, bidegree m, n in the R\, P^. We will show the result for a monomial 
whose degree has increased by one in one variable. Such a monomial can be written X^M 
or MXf^, with M of bidegree m, n (recall that Xj^ — {R\, P^)). We will concentrate on the 
first possibility, the proof being obviously the same for the other case. We compute directly 
X^MiXo, 0) + J^dXiM (X„ a). 

Given our formula for the differential of a product as well as the recurrence hypothesis 
which states that the integral formula is true for M we have : 



X^M(Xo, ^) + j^d (X^M (X„, a)) = " ^ dXl}j (^M{Xl h) - dM (X«, a) 

+ / dX^M (X«, a) + / X'JM (X„, a) 
Jo Jo 

-J2 f dXldXiV^.^M{X^,a) 

In the last line we have chosen (for the sake of simplicity) not to separate the product 

dX^^dX^ into symmetric and antisymmetric part, but to keep the second order terms in a 
compact form. Note again the minus sign arising in front of this term due to our choice of 
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definition for the differential. As a consequence one has : 



X^M(Xo,0)+/ d{Xl^M{Xa,a)) 
Jo 

= XlM{Xl h)- dXi^M{Xl h)-Xi / dM (X«, a) + dX^ I dM (X„, a) 
Jo Jo Jo Jo 

+ dXiM{-Ka,a)+ X^rfM (X«, a) - V / dX^JX^V ^,^M {X^, a) 
Jo Jo j Jo 

Now, due again to the recursion hypothesis, we rewrite A — dX^^M (X^, a) + 
X'^dM (Xq, a) in the following manner : 

A = r dXl\M{Xlh)- rdM{X^,,a') \ + / - f dX^^ dM {X^, a) 

Jo [. Ja \ Jo \ Ja J 

ph ph ph ph ph ph 

= / dXiM{Xlh)+xU dM{Xa,a)- dXl dM{X^,,a')- / dX^dM {X^ 

Jo Jo Jo Ja Jo Ja 

Now, we have to take care about the meaning of the two last integrals. Actually, given 



, one has to 



our definition of the differential, in the first integral dX^^ J^dM {Xa',a') 
consider that a' > a, since our differentials are pointing downward, and as a consequence 
when a' is getting closer to a, the "last" element of integration is dM (Xa+da, 0( + da) — 
M (Xa+da,Oi + da) — M(Xa,a). In the integral dX^^/dM (Xa,a) one has rather to 

consider a' < a. As a consequence, one has : 



J" i a 



dXidM {X^',a') + 



a'<a<h 




dXidM{Xa',a') 



pdXl^ fdMiX^.a)- fdXl 
Jo Jo Jo 



a<a <h 

dM (X^, a) 



the last term in the right cancelling the diagonal contribution that should not appear given 
our considerations just above. Note that this contribution would be negligible for an ordinary 
integral. As a consequence, we obtain : 

A = ^ dXlM{Xi,h)+XijyM{X^,a)- JJ dXldM{X^,,a')- JJ dX'jM {X^, , a') 



a'<a<h 



a<a' <h 



dXlM{Xl,H) + XI j dM{Xa,a)-\^J dX^ j dMiX^.a)-} dX'JM{X^,a) 

Gathering all the terms leads thus ultimately to : 

X^M(Xo,0)+ / d{XlM{X^,a)) = X\^M{Xlh) 
Jo 



+ 



ph ph 

/ dX'JM (X„, a) - ^ / dX^dXiV^i^M (X^,, a) 
Jo A Jo 
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Now, the integral on the diagonal Jq dX^^dM (Xq, a) reduces to 
Jq dXl^dX^V M (Xa,a). Actually as explained before in the definition of the 
differential, the second order terms in dM(X.a,a), multiplied by dX^ contribute to the 

3 

order da 2 which yields a zero contribution while summing over a. We thus end up with : 

X^M(Xo, 0) + /'rf (X^M (X„ a)) = X^M(X^, h) (11) 
Jo 

which proves Eq. (10) ■. 

Now, we interpret Eq. (10) as follows. Starting from the classical regime the integration 
over a leads in continuous manner to the fully quantum regime. This original manner to 
achieve a sort of quantification is very reminiscent of the renormalization group technics in 
which the integration of the high energy modes leads to an effective low energy theory with 
all quantum fiuctuations included. 

Eq. (10) allows us to rewrite formally a function F(X^, ^) depending on the physical 
variables in terms of the same function evaluated at the classical variables (Xq, 0) plus an 
integral of a differential. This introduce apparently a useless complexity since the expansion 
of the right hand side of Eq. (10) around X/j involve a infinite series of the infinitesimal 
increment of the dX-^- By construction, this expansion implies trivially that all terms of 
the series exactly vanishes except obviously the initial term which is equal to the function 
F (X.fi, h) . Although trivial, this decomposition will prove useful later. 

Note again that the integral in Eq. (10) has to be understood as being computed down- 
ward as seen in the definition of (Rq, Pq,) : the starting point is at h and the differentials are 
pointed downward. For example VR„F.(iRa = — VR^F(Xa,Q;) (Ra-da — R^). However, 
this equation has to be taken with some care. Actually the sum over the terms proportional 
to c^Xq being a sum of terms of magnitude it converges only if the sum is discretized, 

and this will be implicitly assumed in this paper. Defining properly the continuous limit is 
out of the scope of this paper. 
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B. Expectation operator 

Having now a differential set up, we aim at defining a linear conditional expectation 
operator £ {.) . We first set the following formula : 

S (Gda) = S (G) da, whatever G (12) 
S {F {Xn, h)) = F (X;,, h) and 8 {F (X,, a)) = F (X;,, h) - T 8 {dF (X,, a)) (13) 



The first equality will allow to define recursively expectations of integrals of a function 
F (Xq,, a) with a deterministic increment da. The second one will make the expectation 
conditional by fixing the starting point X/^ of the path Xq,. The third equation allows to 
define the expectation of an arbitrary function F (X„, a) with the help of its initial value 
F (X/i, h) and the integral of the expectation of a differential whose rule of computation 
is given below. Note that, due to the required linearity of the expectation operator, the 
expectation has to commute with the integration. This property is implied by the third 
equality which yields 8 {dF (X^, a)) = j^8 {dF (X^, a)). By analogy with the stochas- 
tic calculus, Xq, will be seen as a random path whose infinitesimal increments rfX^ will have 
zero expectation. However, again as in stochastic calculus, the expectation of a quadratic 
term like dX^^dX^ can not be set consistently to vanish since it can be written in terms of 
8 {dRi^dPi + dPidRi^) and 8 {dRJPi - dP^dE}^) = -i5'^da, and the last term can not van- 
ish due to the first definition in Eq.(12), and only the symmetric part 8 {dR^^dP^ + dP^dR'^ 
can be chosen to vanish. 

Consequently we set the expectation rules for products at the same "time" a 

^ {i[ dRa}j =^{fl =0' ^{fl (d^^dXjj^. + d^aA^^ = (14) 

and for different times 

(n p \ 
H H dX^AX% =0 for a, ^ a, (15) 
i=l 3=1 / 

We need also the independence of the increments c/Xq, with respect to a function cvahiatcd 
at the corresponding X^ or all the "previous ones" Xq,/, a' >a, which is formulated as 

8 {F {X^,) dZ^) =8{F (X„0) 8 {dZ^) = for a' > a (16) 
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where a — (cci, ccp) is an arbitrary number of values all different and all lower or equal 
to a'. dZot is condensed notation for a product HiLi dZai where the dZ^. can be dX^^ or 
{dX^ldXli^ + dXliJX^l) with ki = 1..6. 

We could impose Eq. (16) directly, but this property is in fact the direct consequence 
of a single condition that we will thus enforce which is the independence of the increments 
dX.a with respect to the initial value X?i of the path : 

S {F (Xa) dZ^) =S{F (Kn)) 8 (rfZ„) = (17) 

That Eq. (16) follows from this single condition is a direct recursive computation. Actu- 
ally, Eq. (16) is trivially checked for F(Xa') a polynomial of degree in (Xq,/,^') since 
£ (ad^a) = (dX^a) = for a constant. If Eq. (16) is true for a polynomial of degree N in 
the variables (X^/, a'), then consider F (X^/, a') to be of degree + 1. As a consequence, 

S (F (X,0 dZ^) = ^ ( (Xft) - £ dF (X^, P)^ dZ^^ 

^ 8{F (Xn) dZ^) - 8 (^£dF (X^, (3) dZ^ 

= -£(^£dF{Xp,P)dZ^^ 
the last equality being a consequence of (17). Now, since 

S {dF (X^, P)) ^ S (d^F (X^, P)dp + J2 ^xf (X;3, /9) dX^ 
\ i=i 

1 ^ 

-i E "^^"^xf (X/J, /3) {dX'pdX^ + dX^dXl) 

all terms in the integral are of degree lower or equal to A", so that the recurrence applies 
(recall that since the integration is downward, in the integral over (3, P > a' > ot and thus 
in the set (/?, a) all elements are different) and S {dF (X^, P)) dZa — 0. Thus Eq. (16) is 
true for polynomials. Since all along the paper we consider functions F(XQ,/,a') that are 
converging series in the variables (X^', a') (seen as classical commuting real variables), the 
density of the polynomial in that space of functions ends the proof. 
A first consequence of the definitions and Eq. (16) is 

S {dF (X«, a)) = £ {daF (X,, a) + {F (X„ a))) da 
= S {DaF (Xq, a)) da 
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This formula allows to compute the integral of E {dF (X^, a)) set previously, as an "ordinary" 
integral since the integration element is only da. 

The main interest of these definitions will show up later on when a function F (Kn, h) will 
be written as an series expansion of iterated integrals along the path Xq,. Indeed, as already 
said before, the expansion of the two terms in the right hand side of Eq. (10) around X?i, 
implies a mechanism of cancellation of all contributions proportional to products of dX^"-.. 
The expectation operator is built to cancel separately the contributions of products of dX^. 
at different times and symmetric products {dX^4X]^. + dX^^.dX^^^ at the same time, coming 
from the two terms of the right hand side of Eq. (10). The rest of the contributions form 
the right hand side of Eq. (10) proportional to da cancel exactly. As a consequence of this, 
an arbitrary function can be written 



V-'o / 

Once again we have apparently introduced some useless complexity, since the function 

F/i (X/j) has been replaced by £^(-F(Xo,0)) which cannot be computed directly, plus the 



duced in both terms and thus ultimately have to cancel . However, this formulation presents 
the advantage of replacing X/j by Xq in the first contribution which will simplify some of 
our computations later. The integration over da will allow to start a recursive expansion in 
the Planck constant. But to do so, we first need to introduce an other operator connecting a 
function cvahiatcd at Xq to the same function evaluated at X^^. This is the role of the EBS 
operation defined below that realizes a kind of connection between the spaces of operators 
at two different values of a. 

A second consequence of our definitions is the factorization of the initial value in the 
expectation : 



This property is usually trivial in probability, but here has to be derived given we have started 
with expectations of infinitesimal increments. It comes from iterating the differentiation of 



Fn (Xn) = 8 {Fn (X;,)) = E fF(Xo, 0) + / dF (X„ a) 





Irrelevant sums of infinitesimal increments have been intro- 



S {F {Xn, h) G (X„, a)) = F (X^) S {G (X„, a)),a<h 
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G (Xq) similar to the chaotic expansion in stochastic calculus. Actually, we can write : 

8 {F (Xn, h) G (X«, a)) = F (X^, h) G (Xn, h) - 8 (^F (X^, h) dG (X«, a) do^ 
given Eq. (16), this reduces to : 

8 {F (Xa, h) G (X„, a)) = F (X^, h) G {Xn, h) - 8 (^F (X^, h) J D^G (X„, a) da^ 

= F {Xn, h) G {Xn, h) - I 8 {F {Xn, h) D^G (X«, a)) da 

J a 

and once again the recurrence works recursively for polynomials. If G (Xq, a) is of degree 
in (Xq, a), then D^G (X„, a) is of degree N — 1, and thus the property can be assumed 

to be true for DaG (X^, a), so that : 

8 {F {Xn, h) G (X„ a)) = F {Xn, h) G {Xn, K) - F {Xn, h) 8 (^J D^G (X„ a) da^ 

= F{Xn,h)8{G{X^,a)) 

as needed. 

We close this paragraph by a general remark about the nature of the expectation oper- 
ation. Despite the formal similarities with the stochastic calculus, the expectation defined 
above is really different, since the expectation is in fact a function of operators evaluated at 
Xn, not a number. This has some consequences concerning the cancellation of certain types 
of expectations that have no counterpart in usual probability theory. But to derive such 
consequences, we first need to define some more tools. 



C. EBS operator 

We now introduce a way to relate two expectations of the same function evaluated at 
two different variables namely 8{F{Xa,a)) to 8 {F {Xa/,a)) but with constant explicit 
dependence in a. Given our previous definitions Eq. (7), one has 

6 

F {Xa-da, oc — da) — F (Xq, a — da) — — {F (Xq,, a — da)) da — x^F (Xq, a — da) dX^^ 

i=l 

1 ^ 

+ 4 E ^xi'^Xi^F (X„, a - da) {dX^dXi + dX^JXQ 

(18) 
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The absence of the d/da term imphes that this expression really depends on the choice of 
the symmetrization of the function F (X^, a). As a consequence considering the expectation, 
one has : 



where (.) is the bracket operator defined previously. At this point we have to remark 
that, since the (.) operation depends on the symmetrization procedure, the right hand 
side of Eq.(19) seems to depend on the choice of symmetrization in the variables in 
F(Xq,,q; — da). However, and this is the most important for us, the operation that sends 
F (X^-da, a — da) to (1 — (.) da) F (Xq,, a — da) by first replacing X„_dQ, by Xq, in F and 
then applying (1 — (.) da), is independent at the first order in da of any choice of sym- 
metrization for the operators in the series expansion of F {^a-da-i — da). This is a trivial 
fact since ((1 — (.) da) F (X^^, a — da)) — F (X^, a) — + (.) da) F (X^, a), and this ex- 
pression is independent of the symmetrization choice as explained before. This operation, 
that will be the only one relevant for us, is thus a consistent and independent of any sym- 
metrization. The apparent trouble with the right hand side in Eq. (19) comes only from 
the fact that changing the order of operators ^a-da is not the same operation as changing 
the order of operators X^ due to different commutation relations. To avoid any confusion 
a subscript a — da should be added to F to remind that the symmetrization has to be 
considered for the variables yia-dai but it would complicate the notations in a useless way 
for the present paper. 

The meaning of our operator l — da being now clarified, the above differential equation 
Eq. (19) can be integrated to get the following relation 
Proposition 



valid for a\> a2. 

The bracket (.)^ reminds that all commutators between canonical variables have to be 
computed with values ia. The operator 5'x„ is a shift operator that sets the dynamical 
variables to X^. that is S^^F (X^) = F (X^). Its action does not depend on the relative 
values of a and (5 and (S'x^-S'xo = S-j^^ whatever the values of a and (5. In addition, T is the 



£ {F (Xa-da, OL - da)) =£{{!- (.) da) F (X^, a - da)) 



(19) 



£(F(Xc,„a2))-^ -Sx„^ Texp 




[^x„,(.)«^xjc^a F(X«„a2) (20) 
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notation for the ordered exponential, that is : 

oo „ 
"-%</3„<../3i<ai 



/ 

"-%</3„<../3i<ai 

Once again, the result is independent of the choice of initial symmetrization chosen in the 
variables X^j. The introduction of the shift (S'xo is crucial and induced by our previous 
considerations since the bracket operation has to be accompanied by a shift of the variable 
all along the integration process, so that ultimately is replaced by Xq,^. 
Proof. 



The proof of Eq. (20) is as follows. Define the function : 



The dependence in the parameter a2 has been skipped for the sake of convenience. Differ- 
entiating G (Xai, tti) with respect to ai yields : 

dG{Xaj^,ai) = G {Xai,a-i) - G {Xaj^^da,^! - da) 



—G (Xo-j-da, «! — da) 
(l-(.)^^dQ;)5x„^5x„^_,^^ 
-G(Xc,i_da,Q;i - da) 



where the last equality is taken at the lowest order in da. Considering the expectations of 
the quantities involved, and using Eq. (19), we are thus ultimately left with : 

EdG (X^^, ai) = ^ ((1 - {■)ai da) Sy^^_^G {X^^_da, «! - da) - G {X^^_da, - da)) 
= ((1 - da) G (Xq^, ai - da) - G (X^i-da, - da)) = 

and as a consequence, using the properties of the expectation operator : 



a2 



8{G{Xa^,ai) - G {Xa^,a2)) ^ S / (X^^, ai) = / ^dC (X^^, ai) = 
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The fact that G (X^j, 0:2) = F (X^j, (^2) ends the proof ■. 

To gain some space in the sequel, we will define an abbreviation for the previous operation. 
We will denote ultimately : 

exp (- (.>1^„,) ^ Texp {- £ ^x., S^Jo^ (21) 

as the exponentiated bracket plus shift, called EBS operation. As a consequence, Eq. 
(20) rewrites as : 

S (F (X„„ as)) = S (exp (- F (X„„ as)) for > ^2 (22) 

This last formula can be understood intuitively as follows, since the EBS operation changes 
the function (by the action of the bracket defined above) but also changes progressively the 
variables from X^j to X^j (through the shift operator) both the bracket operation and the 
shift of variable compensate each over to produce the equality with the left hand side. 

Let us ultimately insist on the fact that the expression in the right hand side 
exp ^— {■)a2^ai^ F (Ka^, is a function of (Xq,^,q;i) due to the action of the shift and 
is therefore a different function than F (X^j, 0:2) . Only the expectations of both expressions 
of Eq. (22) are equal. 

In the particular case where cti = h, since the right hand side of Eq. (22) depends only 
on (X^, h) the expection can be safely removed so that 

£ {F (X„„ 0^2)) = exp (- F (X„„ 0^2) (23) 

We end up this section by giving a full generalization of the previous formula to a product 

of n functions and show the 
Proposition 

S {Hi (X„i, ai) ....Hn (X«„, an)) 
^ S (e-^ )k^'^Hi (X„,, ai) e-W'^^-i/Z^ (X,„ a^) ...e-<->«»— (X,„, a^)) (24) 

for an < ... < a\ < a. This expression for a — h becomes 

8 {Hi (X 

an Q!l) — Hn (Xq,^, an)) 

= S (e-<->«i-'^i/i (X„,, ai) e-<->'2^«ii/2 (X«„ 0^2) ...e-<->--«-ii/„ (X„„, a„)) (25) 
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Note that the order of the functions in the product is irrelevant, and the formula is also 
true if some functions are permuted. 

Proof. We start the proof with the following lemma: 
Lemma 

Let Fi (Xq,^, ai), F2 (Xq,2, ai) be two arbitrary functions with a2 < ai- We state that : 
£ {dF2 (X«„ as)) = ^ £ (Fi (X«,, ai) dF2 (X«„ «2)) = (26) 

This result has no counterpart in usual probability, and, as said in the previous section, 
this is the consequence of our definition of the expectation which sends an operator to an 
operator, not to a number. 

Proof of the Lemma. 

We first show that, ior a < H and every function F (Xq, a) that can be expanded in 
a series of (Xq,,q;) converging "in the classical sense, that is when the (Xa,ci;) are seen as 
commuting real variables, S {F (X^, a)) — implies that F (Xq, a) — 0. This is trivially true 
for a polynomial of degree or 1 in X^ since X^ = X^— J^dX^ and £X„ = X^— J^£^(iX^ = 
X.fi as derived from Eq. (16). As a consequence £ (aX^ + b) = (aX/j + b) which is null only 
if a = 6 = 0. 

Now, assume that S {F (Xq, a)) — F (Xq,, a) — for all polynomial of degree lower 
or equal to N in the variables Xq,. Consider then for F (Xq, a) a polynomial of degree A'" + 1 
(the coefficients depend on a). Then using the EBS operator, we can write : 

£ {F (Xq, a)) = £ (exp (- F (X^, a)) 

given the bracket operation is defined through a derivative of second order, we thus have : 

£ {F (Xq, a)) = F (Kn, a) + Polynomial of degree N - 1 in 

If this is 0, then necessarily the monomials of degree + 1 in F (X^, a) have to be null. As 
a consequence F (X^j, a) is of degree A^ and thus null by hypothesis. As a consequence, the 
lemma to be shown is true for every polynomial, and by a density argument for all the kind 
of series considered. 

In a second step, notice that a direct consequence of the above proposition is that 



£ {dF (Xq, a)) = ^ D^F (Xq, a) = 
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actually, £ {dF (X^^, a)) — £ {D^F (X^, a)) da — 0, which implies that D^F (X^, a) — 0. 

We can then prove the required proposition. Actually, by construction of our expectation 
operator, since 0:2 < ai, 

£ (Fi (Xqi, ai) dF2 (X^^, Oi2)) 
= £ (Fi (X,,, ai) {d^,F2 (X«„ a^) da^ + V x^„/2 (X„„ ^2) dX^^j^ 

+£ [fi (X„,, ai) (^-Iv^.^ Vxi^F2 (X«„ 0^2) dX^^dX^ + dX^dX^^ 
= £ {Fi (Xa^, Oil) -0^2-^2 (Xa2, a2)) rfci;2 
= 

the second equality is a consequence of Eq. (16). This last result shows the lemma ■. 

The proof of the proposition is now in two steps. We first consider a product of two 
functions and show : 

£ {Hi (X,,, ai) H2 (X,„ ^2)) = £ (e-<->'i-i/i (X,,, ai) e-'^-^'^^'^^H^ (X,„ ^2)) 

for ai > a-i (27) 

To do so, let us start again with G (Xq,^, 7) = exp ^— (^)^^^_^^ H2 (Xq,2, a2), 7 is an arbitrary 
parameter varying between ai and 0:2. The dependence in 0:2 in G (Xq,^, 7) is forgotten here 
for the sake of simplicity. By differentiation with respect to 7 we have, as before, at the 
lowest order in da : 

£ {dG (X„^, 7)) = ^ {H2 (X„^, 7) - H2 (X,,_,„ 7 - ^^7)) = 

Now, multiply by Hi (Kan ^1) take the expectation. Since 7 < cci, one has using Eq. 
(26): 

£{Hi{X^^,ai)dG (X,^,7)) -0 

the integration of this last relation (recall that the integral and the expectation commute), 
yields : 

£ (^Hi {Xa„ai) £ dG (X,, a)^ = 

that is : 

£ (Hi (X„,, ai) e-<->'2^"ii/2 (X„„ aa)) = £ {Hi (X,,, ai) H2 (X,„ ^2)) 
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Then, use Eq. (22) with F (X^^ai) = (X„^, ai) exp H2 (X^^, 0:2) between 

«! and h to get : 

£:(//i(X„,,Q;i)//2(X„„a2)) = £:(i/i(X„,,Q;i)e-<->«2^«i//2(X«„Q;2)) 

= £ (e-<->'i-ifi(X„,,ai)e-<-^^2^<^ii72(X,„a2)) 

which is the required result. In a second step, the generahzation to a product 
of n arbitrary functions is shown recursively by starting from the right and re- 
placing i7i(X„„,a„) by exp i7„ (X^^, «„), then to do the same with 
Hi (X„„_i, an-i) exp (- H^ (X^^, «„) and so on ■. 

We can now exploit the mathematical construction developed in this section to consider 
the formal diagonalization of an arbitrary matrix valued quantum Hamiltonian. 



III. THE DIAGONALIZATION PROCEDURE 



We now consider a generic matrix valued quantum Hamiltonian H (R, P) where R and 
P are the usual canonical coordinate and momentum operators satisfying the canonical 
Heisenberg algebra. Our goal is the find an unitary transformation U such that U HU^ is a 
diagonal matrix valued operator e (R, P) (block diagonal for the Dirac Hamiltonian). This 
is in general an excessively difficult mathematical problem. For this reason we consider this 
problem by dividing it in several steps. 



A. The Hamiltonian 

First we introduce the unitary matrix Ua (X^) = U (X^, a) which diagonalizes the Hamil- 
tonian H (Xq) where the canonical variables X^ = X have been replaced by the running 
ones Xq,, so that we can write 

(X,) H (X,) C/+ (X,) = (X,) (28) 

with Ea (Xq) = Ea (Xq, a). With the help of the previous identities Eqs. (10) (12) (14), we 
can compute En (X^) = e (X). Indeed as e (X) — £ (e (X)) we can write 

e{X)^£ (^eo (Xo) + ^ dEa, (XqJ^ . (29) 
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Clearly Eq (Xq) corresponds to the diagonal representation of the original Hamiltonian H{X.) 
where the canonical operators X have been replaced by the classical variables Xq. In practice, 
it is usually quite easy to diagonalize the Hamiltonian when the operators are commuting. 
This is the essence of the method. Starting with classical variables we can recursively 
introduce more and more "quantification" through the running parameter a until we get 
the full quantum Hamiltonian. This procedure is now described in the following. 

The quantity S {dSa-^) — S (^(^-^ + (■)^ ^ai (^aj^ dai can be straightforwardly com- 
puted by using the rule previously given for the bracket of a product Eq. (9). Indeed we 
find 

- \ (Vp,C/„,Vk,//(X,J C/+ - V«,C/„,Vp,//(X«J 

- [Vp^U^.H (X«J Vh^UI - VR^Ua.H (X«J Vp,t/+ ) (30) 

As in [8] we now introduce the notations — zC/q^ (X^J Vp;C/+^ (X^J and — 
—iUai O^ai) ^R^Ua-^ (XaJ SO that Eq. (30) can be written in a more useful form in terms of 
quantities of physical interest. Actually, replacing H (X^J by ^/^^eai^ai everywhere, allows 
to write 

iVp.U^yR.H (X„J UX - Vn.U^yp.H (X„J U^) 
{Ua^P^H (X«J Vp,C/+ - U^^rM (X«J Vp,C/+ ) 

and 

{VpV^.H (X.J VrVX - Vr.U^,H (X„J Vp,t/+ ) 

- -A^'F A^' --A^'F A^' 
~ 2 "1 "1 "1 2 "1 1 "1 
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so that ultimately, one has: 

(31) 

We note that by construction dSai is a diagonal matrix and we obviously have the fol- 
lowing identities 

(a^ + (•)) (X«J = P+ (R.H.S. of Eq. 31) (32) 
= P_ (R.H.S. of Eq. 31) (33) 
where V+ and V- are the projection on the diagonal and off the diagonal respectively. 

+ (•)) (XaJ = (X„J + 7^+ (L>„ (X«„ ai) // (X,J) C/+ } (34) 

with is given by the projected "linear part" of the r.h.s. expression of Eq. (31) : 

Oai£ai (XaJ = V+ ((-Da (X^i, Oi) Uaj_) U'^^Eo,^ + Eo^^Ua^Da (X^i, ttl) t/^J 

Let us remark that the operator 0^^ is well defined since, as explained before, 
the operation Da{y>^ax-,Oii) = (^af^ + (•)) independent from any symmetrization 
scheme for e^i (X^J and if(XaJ, and so is O^^.^aj (X^J = (X^^, ai) e^i (X^J - 

V+ {t/ai (i^a (Xa„ «i) (X„0) C/+ }. 

Since in the applications of practical interest (Bloch electrons, Dirac Hamiltonian...), 
{Da (Xq,j, Oil) H (Xq,J) cancels, we will set this term to for the sake of the exposition and 
will consider later its contribution. As a consequence, one has : 



£ (X) = £ (^£0 (Xo) + ^ Oai.£ai (X^J da^ . 



h 

(35) 
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Now, similarly to Eq. (29) we can write : 

(X,J = So (Xo) + de^, (X,J = So (Xo) + / d (C/ (X,J H (X,J C/+ (X,J) (36) 



which can be inserted in the expectation Eq. (29) to get the full quantum diagonal repre- 
sentation e (X) = £/j (X^) as : 



£h (Xa) = S (^So (Xo) + dSa, (Xc,J da^ 



S So (Xo) + / O (XaJ dai 



s(^So{Xo) + J^ 0«,. So{Xo) + d[U{X^,)H{X^,)U+{X^,)] 



dai ) (37) 



Now remark that since a2 < ai, the terms Vx„ iXa2) dXc 



and 



V^.^Vx.^£„,(X„J(rfX;^rfX^^ + rfX;^rfX^J in de{X^,) = d [[/ (X„J i/ (X,J f/+ (X„J] 
do not recombine with anything coming from Oai (X^J to give a product of the type 
(iRajCiPajthat would induce a ciQ;2 contribution to the expectation. They will thus cancel 
in the expectation. As a consequence, Eq. (37) can be written: 

SniXn) = S ^£o (Xo) + Oa^. So (Xo) + Oa^-Sai (XaJ da2 dai^ 
Repeating the procedure one can then show by iteration that : 



en{Xn)^S 



1 + 



n=l •^0<a„ 



Oai...Oa„dai...dar 



So (Xo) 



(38) 



<...<ai<?i 

At this point, two important comments have to be made, both related to the fact that 
the operators Oa^ have been designed to depend on Xq,.. First, notice that the gradient 
appearing in the definition of Oq. have to be taken with respect to Xq,., but that these 
operators act on functions of X^j with aj < ai, such as £o(Xo). This is not a problem 
however, since Xq,^. = X^. — f^^ dXx and thus the derivative of a function F (X^^) with 
respect to X^^ is just the same as the derivative with respect to X„^. , that is the gradient 
of the function. This the reason why, in the definition of Oa^ we discarded any reference to 
the Qii's in the gradients. 

The second remark is that we can use the EBS operation defined in Eq. (21) and formula 
(25)to put the series appearing in the previous formula Eq. (38) at the same point. Once 
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again, since < q;„ < ... < cci, we can write, inside the expectation : 



Oa....Oa„dai...dan 



0<a„<...<ai<H 



So (Xo) 



p \-/ai^nf) p ^•'«n-l^"n-2 O P 



\ / a.n—*c 



0<an<...<ai<?i 



£o(Xo) (39) 



Note that this formula is a generalization of Eq. (25) since the Oa,. are not functions, but 
rather operators including derivatives with respect to the canonical variables. The part of 
Oa^ acting through some multiplication on the left or on the right or both on the left and on 
the right is not problematic and Eq. (25) applies (with the slight modification of an irrelevant 
change of order in the multiphcation of the functions with respect to Eq. (25)). Likewise, 
the differential part of Oq,^ (Xq^.), such as \A^[ (X^^.) Vr;^ acting on some function say 
F (Xaj._|_i), is not problematic since the operator cxp ^— is a (complicate) func- 

tion of the derivative with respect to the canonical variables and thus commute with Vri 



As a consequence setting \A^l (Xq,^) Vr;^ F {Kak+i) the same point through the EBS 



operation , that is replacing it by \A^[ (X^J exp [- ^ 



in the 



expectation really amounts to compute 



MS (XaJ V«< ^ exp (- F (X„,,J 



that is to let the EBS operation act before the action of 0^^ i^au+i)- This thus justifies 
formula Eq. (39). 

The computation of Eq. (39) derives directly from the definitions of the various operators 
involved and proceeds as follows : for n > , starting from Sq (Xq), first apply e~^'^o^«" on 
Eq (Xq) that will shift Xq by Xq,^. Then make Oq,„ acts on the result, apply e ^ '""^"n-i which 
replaces Xa„ by Xq,„_^ and so on. Ultimately at the end of the process, Xq is replaced by 
X^, this last variable being independent from the integration variables, the various integrals 
can be computed easily, as the integral of a series expansion. This process of integration 
will be of course the same for all similar expressions in the sequel. 

The case n — Q corresponding to the first term in the series expansion : S^Sq (Xq)) is of 
course easily handled by replacing it with e~^-^o-^'i£o (Xq). 

We can now deduce that (with the notation X = X^) : 



£(x) = £n(Xn) = £: 



Texp 



£o(Xo) 



(40) 



.JO<a<h 

which is a compact expression for the diagonal Hamiltonian e (R, P) in terms of the "clas- 
sical" diagonal Hamiltonian Sq (R, P) in which the classical variables Rq, Pq have been now 
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replaced by the quantum ones R, P due to the EBS action. Here T is the usual notation 
for the "time ordered product". Eq. (40) is the required expression and constitutes the 
main result of this paper. 

Wc can now consider the additional contributions that would appear if we 
had considered a case such that (Dq, (Xq,^, ai) if (X^J) 7^ 0. Let C (X^ J = 
P+ {[/q;i (Da (Xq^, cti) if (Xq,J) In such a case, the repeated application of + 

C (Xq,J yields rather : 



en{yin) = £ 



1 + E 



Oai---Oa„dai...dan 



^—l J 0<a„<...<ai<a 



So (Xo) 



1 + 



00 p 

u 



Oai ■ ■ ■ Oa„-i dai.. .dan 



C(X„J 



' 0<a„<...<ai<a 

Using the same tricks as before it leads directly to the following expression for the diagonal- 
ized energy operator : 



+S 



Texp 

h 



ceo 



0<a<h 

Texp 



£o(Xo) 

— / \^ —( \^ 



ao<a<H 



C (Xq,J dao 



Clearly the practical apphcation of Eq. (40) requires the knowledge of the transformation 
matrices Ua which enter into the definition of the operators Oa- 



B. The transformation matrix U 



Note first that there is a certain arbitrariness in the choice of the unitary matrix Ua (X) 
as explained in [8] which reflects a kind of gauge invariance. Actually, multiplying the 
transformation matrix Ua (X) on the right by a diagonal unitary matrix yields an other 
diagonalization, equivalent to the previous one. In particular it allows to choose n conditions 
for the diagonal entries of [/„ (X) {n x n being the size of Ua (X)). An exphcit choice will 
be done below to simplify our expressions. 

To find the transformation matrix U we use the same approach as for the diagonalization 
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of the Hamiltonian, by writing : 



Uh (X,) = £ (u (X„) + £ dU„ (Xj") 



we can find {da + {■))Ua (X^) by using again the diagonahzation process leading to Eq. 
(31). Indeed we had 

(42) 

In addition, the unitarity condition UaU^ = 1 implies the relation {da + (•)) {UaUa ) = 
which reads : 

= ((1^ + (•)) u}j U+ + Ua(J^ + (.)) - I {Vp^UaVn^K - Vn^UaVp^U^) 
or, in a more compact way: 

Mixing the two equations Eqs. (42) (43), we obtain that Eq. (42) after projection on the 
non diagonal part becomes the equality 



V- 



- I'P- {^5 Vfi,£„, + VR,ea,A^[ + A^[Vp,ea, + Vp^s^.A^ } 

- -V iA'^'E A^' -A^'£ A^' +£ \A^' A^'^ + \A^' A^' 

(44) 

Both conditions Eqs. (43) (44) can be solved for {da + {■))Ua- Actually, decomposing 
[{da + (•)) Ua] Ua = hr + tthr in Hermitian hr and anti-Hermitian ahr part, Eq. (44) reads 

V- [hr + ahr, = M {Ua) (45) 
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where M (Ua) is the r.h.s of Eq. (44). We can now use a particular choice of gauge. We fix 
our diagonahzation process by setting V+ (ahr) — 0. Practically it corresponds to multiply 
Ua on the right by a unitary diagonal matrix Da (which has n degree of freedom) such that 
V+ ((^ + (•)) UaDa) D^U^ is Hermitian. This gives n conditions that determine Da- 

Now, from Eq. (45) and given our gauge choice for Ua, we deduce that V- {ahr) = ahr 
and [ahr, Sa] — \ {M {Ua) + M^{Ua))- To find ahr we aim now at inverting the commutator 
[., £q,] in the left hand side of this equation. First note that, since we have assumed the diag- 
onahzation is possible, we can safely assume that the right hand side | (M {Ua) + M+(C4)) 
lies in the image of the operator [.,£„]• It is in fact the case in all practical cases. Since both 
ahr and M {Ua) are anti-diagonal by construction, we just need to invert the commutator 
with Sa, on the space of anti-diagonal matrices. To do so, we first need to find the 

kernel of [., Sq] on this space and check if it is null. 

Assume first that the bands of the system are non degenerate, that is (£«)„ 7^ {^(y)m 
m^n. As a consequence an element V in the kernel of the commutator satisfies 

Km [(£a)^ - (£a) J = (Vp, (£«)„ Vp^Km " Vr, (£«)„ Vp.Km) 

that is : 

Km = (Vp, (£a)n Vfl^Km - Vp, (£«)„ Vp.Km) r, . \ , . . 

Since the beginning, we have assumed that all the functions involved are regular in a and 
have a series expansion in this parameter (as well as in the canonical variables). Thus, 
iterating the last relation for Vnm yields that, for each Vnm expanded as a series expansion 
in a, at each order in a, Vnm = 0. As a consequence, the kernel of the commutator with Sa 
is null and we can thus formally write ahr under the form : 

ahr = TV (a) .Ua (46) 

where we have defined : 

^ («) -Va = \ ([•,£«]"') (M {Ua) + M+{Ua)) (47) 

and the inverse of the commutator operation [, satisfies obviously : 

[[., £q;]~^ .M, £a\ — [.,£cX'^ ■ [M, £q] — M for M anti-diagonal 
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If the Band are degenerate, we need to define tfie inverse of the commutator more carefully. 
Actually since some bands may be degenerate, the kernel of the commutator [.,£a] has no 
reason to be null now. This operator is thus no more bijective and we will not be able to 
define its inverse uniquely. This non uniqueness will lead to a non unique definition of ahr. 
This is not astonishing. Actually, when the Band are degenerate, as in the Dirac case, we do 
not look for a diagonal Hamiltonian,but rather for a Block diagonal Hamiltonian. This leaves 
thus more freedom for the diagonalization process and as a consequence, for the definition of 
[.,£a]~^- Practically, we proceed as follows. Decompose our space of non diagonal matrices 
of operators denoted AD in the following way. Write AD = ker © adi with adi an 

arbitrary supplementary space of ker ([., Now, given M an clement of Im ([., Sq]), chose 
a antecedent mi of M with respect to [.,£«], and decompose uii with respect to the previous 
decomposition, mi = 1112 + 1713 where m2 € ker ([., ej), e adi. Now set [., ej"^ — ^3- 
Given a chosen decomposition ker ([., ea\)®adi, is unique since if mi is another antecedent 
of M, rhi — mi e ker([.,£Q;]) and as a consequence rhi — rris + (m2 + mi — mi) is the 
decomposition of ttii. 
One can check that 

= [mi -m2,€a\ 
= [mi,£a] = M 

as needed to solve Eq. (45). On the other hand, notice that : [.,£0]"^ [M, is not equal to 
M in general. Actually, choosing M as an antecedent of [M, Sa] and decomposing it gives 
an element rris that a priori depends on the supplementary space chosen. As a consequence, 
the operation [.,£cj]~''^ is only a right inverse to [.,£a]- This is not a problem, since it allows 
anyway to write our solution for Eq. (45) as : 

ahr = N (a) .Ua (48) 

where we have again defined : 

N {a) .U^ = i {[.,6^]-') (M {U^) + M+{Ua)) (49) 

The non unicity in the definition of [.,£a]~^ is reflected in the choice of decomposition for 
the space AD. Let us insist on the fact that this choice is itself the consequence of the 
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Bands degeneracy that allows for a larger freedom of gauge choice for a block diagonal 
representation. 

To complete the determination of Ua we still need to deduce hr which is readily obtained 
from the unitarity condition Eq. (43) as: 

hr^-\[A^^,A^^] (50) 

Indeed, Eq. (50) follows easily from Eq. (43), when noting that ((^ + (.)) Ua) U+ + 
Ua + (■)) is hermitian due to the fact that for any operator A, (.) = ((.) A)'^ (and 
this is true in particular for U^)- This assertion can be checked on an arbitrary expansion in 
R and P of A, since for an arbitrary monomial RiPj'An^k with ^ an arbitrary coefficient 

matrix, (•) (i^F^z'^Xfe = hi^kR^'Pt') K,k = (-f^^^'-'^r'^n,.)^ and this last 
quantity is equal to [(.) [Pl'Rf) An,k]^ . 

Gathering the results Eqs.(47) (50) allows us to introduce the operator Na operating on 
Ua in the following manner Na-Ua — {hr + ahr) Ua, so that we can write 

d 



with Na given explicitly by the following expression 



Na.Ua = I ([.,£a]"') (M (t/„) + M+ ([/«)) - ^ [A^' , A^'] Ua 

As for the energy diagonalization, this expression can be rewritten in terms of physical 
quantities. Recall that M (Ua) is the r.h.s of Eq. (44), so that 

^{M{Ua) + M+{Ua)) = -lV-{A^[VR,Sa, + VR,ea,A^{+A^\Vp,ea,+Vp,ea,Aa[} 

{ A^'] A^' - A^q A^' + H.C.} (51) 

where H.C. stands for the Hermitian conjugate. Note also that Ua ((^ + (■)) H (Xq,)) t/+ 
is hermitian, as a consequence both of the hermiticity of H (Xq,) and the fact that 
((.)// (X„))+ = (.)//(Xq)+ = {.)H{Xa) (see the assertion below Eq. (50)). Ultimately 
we are led to the following expression for Na-Ua- 
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N U = - [ £ 1""^ 



^ + {.)]h{X^)] uA +!-V-{[ea,A^q A:^ - [ea,A^q A^^} + H.C. 



da ^ 
-IV.{[A2\A'J]U^} 

where Sa O^a) is computed recursively as explained before. 
As for e (X), we can therefore write for U (X) = Ufi (X^) : 



(52) 



U(X)=S 



Texp 



0<a<h 



Uo (Xo; 



(53) 



This expression has a very similar structure as the solution for the energy Eq. (40) except 
that the operator Oa has to be replaced by A^^. 

To end up this section, and given the solution derived for Ua, we can rewrite the operator 
Oq in an simpler form. Indeed starting again from Eq. (31) using the unitarity condition 
for Ua yields : 



d 
da I 



+ (■) (X 



d 



da I 



d 
da I 



oil 



1 



(54) 

Now, recall our gauge condition which states that P+ ((^ + (.)) Ua) Ua is Hermitian, so 



that n 



Uav^o'i antihermitian, as the commutator of two hermitian 
quantities. As a consequence, V+ (^(^^fl ^■'^) ^"i) H.C. — 0. Then, use the 

fact that Ea (Xq) is assumed to be an Hermitian operator, so that + (Xq,J is 

equal to half its sum with its hermitic conjugate. A computation similar to Eq. (51), leads 
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us directly to the expression : 



+ 



P+ { [e^, A^^] A^^ - [sa, A^q ' } + H.C. 
d 



da I 



+ (.)^if(X„J[/+ 



The action of Oa is obtained by skipping Ua^ (^-^ + (•) j H (X^J U^^ from the last term, so 
that : 



OaEa (X„) = -V+ {A^'^R^Sa + "^R^SaA^' + A^'V P^E^ + Vp,£„^a'} 



+ 



P+ { A^'] a:' - ^^'j A^^} + H.C. 



(55) 



Let us ultimately recall, that in most of the physical applications of interest for us (in particu- 
lar for Dirac and Bloch electrons) the term + (.)) if (X^) will cancel. This simplification 
will be assumed for the rest of the paper. 



C. The full system H and U 



We can now write the solution of our diagonalization procedure for a general matrix 
valued Hamiltonian through an unitary transformation U as the solution of the following 
system of differential equations 



£(X) 

t/(X) 



£ 
£ 



( r exp / e"^->«-fiOae"<->o-«rfa £0 (Xq) ) 



Texp 
Texp 



0<a<H 



£0 (Xo 
Uo (Xo) 



(56) 
(57) 



where and are given respectively by Eqs.(55) and (52). The only prc-rcquircmcnt is 
that the diagonal form at o; = 0, Eq (Rq, Pq) is known, i.e., when R and P arc considered 
as classical commuting variables (which means that C/q (Xq) is known). Of course these 
equations do not allow to find directly e (X) , U (X) since those quantities are involved on 
the R.H.S. of these relations. However, they allow us to produce solutions for e (X) and 
U (X) recursively in a series expansion in h. Moreover, and as needed, the results of our 
previous sections show that the matrices e (X) and U (X) obtained through this process are 
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independent of any choice of symmetrization. A particular choice in the way of arranging 
the variables for U (X) = will lead to differently symmetrized, but identical, operators. 
As it will appear clearly later on, having both e (X) and U (X) at order n in and 

reinserting in the exponential of Eqs. (56) (57) allows us to find e (X) and U (X) at order 
n + 1 in ^. But before solving recursively the set of equations Eqs. (56) (57), we first compare 
the present approach with the one developed in article [8]. 

IV. LINK WITH THE DIFFERENTIAL EQUATION OF REF. [8] 

In [8] we developed a different less general approach which led to the differential equation 
(X„ a) = [dJJ^ (X,) C/+ (X,) , (X,)] + i [A^^Vr^s^ + Vr^s^A^^ + A^^Vp.e^ + V 

+ ^a{H (X,)) C/+ - [BaSa - SaBa] - (Sa) } (58) 

(in the right hand side we have skipped the explicit dependence in a for the sake of simplicity) 
which was coupled to the evolution of the transformation matrix Ua (X^) as a function of 
a [8] : 

= 9„C/„(X„)C/+(XJ+C/„(X„)9«C/+(X„) -^iBa- B+) + [A^^A^'] (59) 

where -|Bq, = {Asym \^RiVpJJa (X^)]) C/+ (X^) = (C/q,) C/+.With these two equations Eqs. 
(58) and (59) at hand, the diagonahzation process can be performed. Actually, since all 
quantities are matrix valued and since (X^,) is by definition a diagonal matrix, we can 
separate the energy equation Eq. (58) in a diagonal and a off-diagonal part such that we 
are led to the following two equations 

d 

—Ea (X„) = P+[R.H.S. of Eq. 58] (60) 
oa 

O^V-[R.ll.S. of Eq. 58] (61) 

In [8] it was claimed that those three Eqs. (59) (60) (61) allow us to determine recursively in 
powers of a the energy of the quantum system in question. Actually, the integration over a 
of Eq. (60) gives Sa (Xq) at order n in a when knowing all quantities at order n — 1. By the 
same token, Eqs. (61) and (59) (whose meaning is that Ua (Xq,) is unitary at each order in 
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a) involve daUa (X^), and allow to recover Ua O^a) order n by integration over a. As a 
consequence, the diagonalization process is perfectly controlled order by order in the series 
expansion in a. In [8] we also provided two physical examples at the order h^. 

Now we want to show that our solution Eq. (40) satisfies the differential equation Eq. 
(58). 

To do so, we first rewrite (58) and (59) with the notation of this paper. First, (59), is in 
fact (43). Actually, -| (B« - B+) = (C/„) U+ + (C/+). 

Second, (61) is identical to (44). It can be seen by starting with (61), and noting that in 

(61) 

[dJJo, (X,) U+ (X,) , (X,)] - [B^e^ - sM 

= [daUa (X,) U+ (X„) , (X,)] + (Ua) U+Sa - 6^ {Ua) U+ 

which, reinserted in (61) and projected on the non diagonal subspace through leads to 
(44) after some rearrangements. Note that the term proportionnal to (X^) is missing 
in (61) compared to (44), since it was assumed to be nul in [8]. 

Now, we show that £q, (Xq) as defined in (56), satisfies (60). We will use the hermiticity 
of £a (Xo, a) as well as the previous gauge condition stating that V+ [[^ + (.)) Ua) U^ is 
hermitian. As proved in the previous section, it implies that [((^+ (•)) Ua) U^,ea] is 
antihermitian. We also assume (i/(X„)) = 0, to be consistent with the previous sections, 
but including it would not harm. We are led to: 

■^EaCKa) = ^P+[R.H.S. ofEq. 58 + H.C.] 
oa 2 

+ \ { [Sa, ^f'] - [Sa. A^a' ' [^a, ^'j] + H.C.] - (s.) 

— OaSa (Xq) — {Sa) 



Now, from Eq. (56) we can write the partial differential with respect to h 
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Texp 
= -(.)£n(X,) 



0<ai<h 



^0 (Xo 



\ " 0<a„<...<a2<h 
= -(.)£r(Xr) + 



-\-/a2^hf) p ^•'a„_l-»a„_2/l p \-'a„^a, 



"-'Oa^e <->o-""£o (Xo) 



\ 0<a„<..<a2<h / 
= -(.)£a(X;,) + 0^£^(X^) (62) 

The first equality has been obtained by using that by construction, -^e ^''^i^fi = 
— >S'xfi {■)fiSxn s-iid by remarking that in computation of the partial derivative the canon- 
ical variables and the expectation operator remain unchanged (that is X^j is seen as a con- 
stant). We have given previously the expression for OnSnC^h)- Using moreover Eq. (59) 
one gets directly that {e (R, P) — Eq (R, P)) is a solution of the differential (58). 



V. DYNAMICAL OPERATORS AND COMMUTATION ALGEBRA 

In this section we will see that new non-commuting position and momentum operators 
which have contributions from Berry connections emerge during the diagonalization and are 
more suitably to correspond to physical operators (for the physical discussion of this point 
see [9] [11] [8] [22]). This discussion may be skipped, and the reader can directly move to the 
next section dealing with the general expressions for the diagonalized Hamiltonian to the 
first and second order. 

From Eq.(31) one sees that the operator Oa = {da + (•)) acting on Eq (X^) can be de- 
composed as a sum of a "translation" operator T and a "magnetization" M operator (this 
terminology is explained in [8]) 



O^Eo (X,) = (r„ + M,)£o (xj. 



(63) 
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where the "magnetization" operator acts as 



M^so (X J = { [e^, A2^] A^^ - [sa, A^'] A^^ } + V, 



(64) 



and 



r„£o (X„) = -V+ {A^'VR,e^ (X,) + Vr^SoA^ (X,) + Al^Vp.e^ (X,) + Vp.e^Al^ (X,)} 

(65) 

To inspect the action of , let us skip the magnetization contribution Mq, for the sake of 
clarity, and consider the following relevant contribution for the computation of the diago- 
nalized Hamiltonian : 



£Texp 



_/ \S _/ \S 



£o (Xo) 



'0<a<h 

When developed in series, as explained in the previous sections, it is given by 



n r p 
n .^0- 



^ > / e 

L'/0<an<-<Ql<?i 



_/ \S _( \S _/ \S -s 



£o(Xo) 



(66) 

Recall that each operation of the kind e ^ '"'^"i-i is going along with a change of variable 



Wefind£^Texp J^^^^f^e ^'^"^'^TqC ^-^o-^ado; £o (X) by the following ansatz. We assume 
that £:rexp [/o^^^^e-«^-ftr„e-<-)o-.«da] Eq (X) = So (Xn+Af'^ where Af = (X;^,;^) 
has to be determined and a particular (arbitrary) choice of symmetrization of the variables 
X has been made. Due to the form of the left hand side in the previous relation Eq. (66), 
A'^ is of order h. We choose to write Sq and A'^ such that the powers of R and P are 
in a completely symmetrized form, that is we sum over all equally weighted permutations 
of the canonical variables in the series expansion of Eq (X) (and A-^) . Then, Eq (X + 
is obtained by replacing X by X + A'^ in the series expansion of Eq. Note that, for later 
purposes (•)^£o (X + is of order fi^. To differentiate the right hand side of Eq. (66) 
with respect to h we proceed in the following way. 

On one hand, assuming that Eq does not depends explicitly on h, ^Eq (X + ^4^) is equal 
to : 



dAf 
dh 



Vx£o (X + ^^) 
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Note that in the last expression, the product Vx^o (X + A-^j has to be understood 
in the sense that -g^ replaces X + ^4^ at each place the gradient is acting on the series 

expansion of Eq- 

On the other hand, the same infinitesimal variation can computed through the left hand 
side of Eq. (66) : 

d 



d 



dh 



n=0 



_/\S _/\S 



D<an<---<oii<h 



So (Xo) 



n=0 



/p (■)ai->h'T p ^■^«n-l^«n-2 7~' P ^'^""-'"n-l 7^ p ^)o->q 



D<an<---<Q;i<?i 



D<Q;n<...<a2<?l 

As a consequence, one is left ultimately with : 



£o (Xo) 



£o (Xo) 



[- (.), + Tn\ So (X + A^) = ^ Vx£o (X + Af) 



(67) 



Now consider V+ [A^] and P+ [-^ajj with the following definitions for the "non-projected" 
Berry connections — i [Ua'VpU^] and — —i [Ua'VnU^]. Define also for convenience 
V+ [A^] = {V+ [A^] ,V+ [Al]). Our last equation Eq. (67) is thus equivalent to : 

[A^] . [Vx£o (X + A^) + Vx {AfY (X) Vx,£o (X + ^1^)1 + H.C. - (.),£o (X + A^) 



dAl 
dh 



VxeoO (X + Af) 



where H.C. stands for the Hermitian conjugate. The term Vx^^o (X + A^^ stems for the 
gradient of Eq (X) evaluated at (X + Af). As before Vx {Af)' (X) Vx.^o (X + Af) has to 
be understood in the sense that Vx {AfY repl aces X + A-^ at each place the gradient is 
acting on the series expansion of Eq. To be able to compare both sides of this equation, one 
has to symmetrize both expressions in the same way. As a consequence the left hand side 
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has to be rewritten in the same symmetrized form as the right hand side : 
\v+ [A^] . [Vx^o (X + A^) + Vx {A^Y (X) Vx.£o (X + A^)] + H.C. 



Vx£o (X + Af) + Vx {A^y (X) Vx,£o {X + A^)\j^ + H.C. + C (X) 

This formula requires some explanation. In the left hand side, the multiplication by V+ [-4.^ ] 
is performed half on the left and half on the right, as implied by definition of the transla- 
tion operator. On the right hand side, all expressions are seen as symmetrized in a way 
that V+ \A^\ has been inserted at each place where the Vx is acting (exactly as for 
^^%h^^ Vx£o (X-|- A^))- The {}g is there to recall this full symmetrization in the vari- 
ables. The term C (X) is the correction due to this change of symmetrization while moving 
the V+ [-4.^] inside the series expansion of Eq. By construction, it involves the powers of 
gradients of V+ [^q] and Sq and is of order ti?. In practice, this term depends specifically 
on the problem at hand (that is on the form of the Hamiltonian) and can be computed order 
by order in h. Since later on, we will consider only the order t^, it will turn out that this 
term will be negligible due to an integration. 
We can now write the differential equation : 

^""Wxso (X + A^)^\ {V^ [Al] . [Vxeo (X + A^) + Vx {Aff (X) Vx^so (X + ^^)] }^ 



dh 



+ H.C. + C{X)-{.)^eo {X + Af) 



where all expressions are now symmetrized in the same way. Now, remark that, given our 
previous remarks, C (X) — (.);j£o (X -|- ^4^) is of order fi^. It imphes that, after integration, 
it will contribute only to the third order in h to A-^. As a consequence, neglecting this term 
in first approximation (that will be indeed the case in our applications) we thus deduce that 
A'^ satisfies the following differential equation : 



^+ [^f] + I [^f] -^x) Af + H.C. 



dh 



(68) 



Given that for h — 0, one has Af^ — 0, the solution of this equation can be written recursively 
as: 



0<a<h 



'PA^a] + l{{rA'^a]-^^)Al + H.c.) 



da 



(69) 



/ 5x,n [A^] da+ f 5x4 iK / ^+ [^"i] 

Jo<a<h Jo<a<h ^11 JO<ai<a 



+ H.C. 



+ 



daida 
(70) 
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Recall that y^x^^ is the shift of variable — > X^j = X, so that in the integrals, the variables 
X are inert, only the explicit dependence in a is integrated on. Let us also stress that in the 
previous integrals, the variables Xq, in and have been replaced by X and are thus 
constant with respect to the a integrations. The reason is that in the differential equation 
Eq. (68), X is seen as constant, only h is running. Actually, the equation involves only the 
partial derivative Having now the solution for as a function of the Berry phases 

we can write our solution for the exponentiated action of the translation operator. Dividing 
in two components with respect to R and P, = (>1r, A-p), one has : 



£^Texp 



.JO<a<h 



£o(X)=£o(x) (71) 



with X = (r, p) and r and p are new coordinate and momentum operators corrected by Berry 
connections terms in the following way : 

r = R+^R 

p = P+^p (72) 

justifying the name translation operator for T^. The inclusion of the corrections due to 
C (X) can be performed in the following way. Shifting A-^ by a correction A-^+SA^ with 
A-^ the solution previously found for Eq. (68), gives the following equation for SA-^ : 

i-p+ [Af] . [Vx {SAfY (X) Vx,eo (X + Af + SAf)] + H.C. + C (X) - (.),£„ (X + Af) 
= ^Vxe.(X..f...f) 

Note that C (X) is computed using the derivatives of A^+5A^. Expanding SA^ (X) as an 
h series expansion of completely symmetrized function SA'^^^'' of R and P (starting with 
n = 3) allows, at least theoretically to find the SA'^^^^ recursively by solving an equation of 
the kind : 

^ F (X, 8Af-^ (X) , 5Af-^ (X)) 

where F (^X, 5^^^"^ (X) , 5^^^"^ (X)) is determined by replacing eg, A^, V+ [A^] by 
their expressions at the n-th order. As said before, in our practical considerations the term 
5A^ (X) will always be neglected. 

The variables x = X + ^4^ we have defined seem thus to be the natural ones arising in 
the diagonalization process. 



41 



Note that at the lowest order we have : 



r = R + hV+ [Af] = R + 
p = P + nP+[<] =P + < 



(73) 



where and Aq are the usual "semiclassical Berry connections" defined previously [8]. 
From equations Eq. (72) we readily deduce the following non trivial algebra 



[p,,pj] = zh'e^ = -if? {VrAp, - VrAp:) + [Ap^, ^pj 

[pi, Tj] = -ihSij + i^'ef; = -iMij - ih^ {VR^AR. + Vp.Ap^) + [^p., Ar^] (74) 



where the terms ©jj are the definitions of Berry curvatures. Of course these non trivial 
commutation relations also give new contributions to the equations of motion and thus lead 
to new phenomena [8] [11][9]. The commutation relations are valid to any order in h, but 
in practice we can compute them as well as the energy e (X) in a series expansion in h. 
Relations Eqs. (64) (71) will be helpful when writing the exphcit expression of e (X) in a 
series expansion in h in the following section. 



The exact expression Eq. (40), can now be expanded in a scries expansion in h. Note 
that we will always identify £/j(Rft, P/j) with e(R, P). Wc also implicitly assume for con- 
venience that all expressions are symmetrized in R and P, in such a way that for all 
expression depending on R and P, in the series expansion of this expression, all pow- 
ers of the momentum are put half on the left and half on the right. Any other choice 
of symmetrization would be of course suitable. Recall at this point some previous no- 



tations : A^ = 2[f/,(X,) Vp„f/+(X„)] and A^ = [f/„ (X„) VR„f/+ (X,)].as well as 
A^ = {A^,A^) and .4.^ = Af. We will also denote Af", and Aq" the zeroth order 
Berry connections evaluated at X^, that is A^" = A^" (X^) = i [Uq (Xq,) Vp^Uq (X^)], 



VI. SERIES EXPANSION IN h 




(XJ = -i [Uo (X„) Vp„C/o+ (XJ] , and A^" = (^^, ^^"). 



A. 



First order in h 



At the first order we obviously get the following expression 
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e (X) = £o (X) + [Oa£o (X,)] da 

Jo 

This expression requires some explanations. The action of Oa on sq (X„) leads to an ex- 
pression F{X.a, Oi) which depends both on Xq, and a. The shift operation replaces X^ by 
X^ = X so that the integration is (trivially) performed on the variable a only. To stress 
this point, let us first introduce a convenient notation. Recall that at the lowest order : 

O^so (X,) V«,£o (X«) + V^,£o (Xc.) A^^ (X,) + A^^Vp^So (X,) + Vp,£o (X«) A^^) 

P+ I [so (X.) , A^q A^^ - \ [so (X.) , A^^] A^^ + H.C. 

where the Berry phases can be replaced by their value at zeroth order (evaluated at Xq,), 
that is : 

0„£o (X„) = P+ U (A^-VroiSo (X„) + Vr„,£„ (X„) AS°' + X"' Vf,£o (X„) + Vf,,£o (X„) A*"' 



^0 (Xa) , Aq"^ 



An 



£o(Xq),^^- A^-'+H.C. 



since [f/„ ((£ + (.)) H (X,)) [/+] = 0, at that order. 

Now, let define the "covariant derivative" operator :T>^^ = (I^r„,I^p„) 



and as usual "Dx = ^x^, so that we can write directly : 

OaSo (X«) = [(Px„£o (X,)) ^^ + H.C] 
Therefore, at this order we have 



£(X)=5o(X) + ^ ^x. lv+[{V^^eo{X^))A^-+H.C.] 

= So (X) + 11^^^+ [(^x£o (X)) + H.C] da + 0{h^) 
= £o (X) + ^7^+ [(Px£o (X)) Al + i/.a] da + ©(n^) 
= £o (X) + [(Px£o (X)) Al + i/.C] + 0(^2) 



da + 0(/i^ 



(75) 
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In the three last hnes, we have introduced a shght abuse of notation. Consistently with 
our conventions, we should have written j^'^ rather than A^. However, since in the final 
expression, all terms have to be expressed at X^^ = X and not at Xq, no confusion should 
arise and will, from now, always stand for A^'" in the final results. Note that, as stated 
before, the integral over a in the previous computations is trivial once the shift of variable 
has been performed. 

The previous formula can be expanded in several ways, one has. 



e (X) = So (X) +r+i^ [a^'Vx.So (X) + Vx^^o (X) A^') 



£„(X),.4?'U„'''- £o(X),.4; 



(76) 



or, in terms of covariant variables: 



em 



eo (x) + -V, 



So (X),<' 



A, 



Pi 



£o (X) , A, 



An 



One then recover the formula first derived in ref. [8] where x — (r, p) are given by the 
expression Eq. (72). As previously mentioned, we only need to know Uq (X) and thus the 
zeroth order Berry phases to get the energy expansion at the first order. 

We end up this paragraph by providing the first order expansion for U that will be used 
for the second order diagonalization: 

U (X) = C/o (X) + r S^, [NM (X«)] da = Uo (X) + hU, (X) Uo (X) 

^0 



= (l-h[.,eo]-\ 



r-i\ [a^^Vr^s^ (X) + Vij,£o (X) A^^ + A^,^Vp,eo (X) + Vp,£o (X) 



{ [e, (X) , A^] A^' - [so (X) , A^,] } 



aRi An 

•^0 ' -^0 



t/o(X) 



(77) 



B. Second order 



Given our choice of symmetrization {sq (X^^))^ is of order 1 in a, so that. (■)o-^.a -^o (X) = 
/o'^'S'xa (^0 (Xa))Q,rfQ; is of order 2. As a consequence, expanding the compact form of the 
energy operator to the second order, leads to : 



(X) = So (X) + / S^, [OaSo (X,)] da+ / S^, [Oa,S^^, [Oa.£o (X„J]] da2da, 
Jo Jo Jo 

- / -^x, {eo (Xa))^ da 
Jo 
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The last contribution can be computed easily. At the lowest order in a, ^Sx^ (sq (Xq))^ = 

I {sq (X))^ (actually, S^x.^ (^o 0^a))a ^ function of X times a and (sq is the same 

function times H). As a consequence, /q^S'x^ (^o (Xq))^ (io; = jj^ ^ (^o ~ f (^o (X)) 
(with the convention (ero (X)) = (^o (X))^). 

The first contribution Oadaso (X^) can be expanded as before as: 

/ 'S'x^ [Oa£o (Xa)] da 
Jo 

= ^ 5x,P+ { I (^a' Vk,£o (XJ + Vr^So (X,) ' + Vp,£o (X«) + Vp,£o (X„) 1 da 
+ fsx,V^ [\ [e, (X.) , A^^] A^^ - \ [so (X.) , A^] A^^ + H.C. 



+ 



9a 



+ (.) i/(x„) U/, 



4 
da 



but now, the Berry connections have to be expanded to the first order in a. Note that in 
the second order contribution, due to the double integral, 

J J 'S'x;, [0«i5'x„^ [Oa^^o (Xaa)]] rfas^ai - - (£o (X)) 

the Berry connections need only to be expanded to the zeroth order in a. 

To go further, we thus need to give first some expanded formula at the first order in a 
for the Berry phases A'^ intervening in the definition of the diagonalized Hamiltonian and 
the dynamical variables. 



1. Berry connections at the first order 

To obtain A'^ we only need the matrix U at the first order as computed before : U (X) 
Un (X) = Uo (X) + hUi (X) Uo (X) with Ui (X) given by (Eq. (77)) 



U^{X) = -[.,eo]-'. 

+^{[^o(x),A''' 



P_|i(A^'VH,^o(X) + V 



R^eo (X) + A^'Vp.so (X) + Vp,£o (X) A^' 



An 



So (X) , A, 



A^' 



A^' A^' 



(78) 



At the same order the (non-diagonal) Berry connections A'^ = {A^, A^) are again given by 

A^ (X,) = i [U^ (X,) VpC/+ (X„)] = (X„) R«C/+ (X„) 

and Al (X«) = -i [Ua (Xc,) VR.C/+ (X,)] = -U^ (X,) P«[/+ (X„) 

a 
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where now O^a) is the transformation to the first order in a i.e. Uq (X^) + aUi (Xq) , in 
which X is replaced by the running operator X^. Using the hermiticity of A^, so that one 
has = (^A^ + {A^y^ /2 we can expand A^ as : 

= Q ([1 + aU, (X«)] Uo (X«)) ^ (C/o+ (X.) [l + aU^ (X„)]) + ^H.c}j - ^ 

(the ^ factor reminds that in our definition of A'^ the gradient with respect to X^^ is 
normahzed, i.e. divided by a). After some recombinations, the previous expression can be 
written in a more convenient form : 

= ^Uo (X„) [X^,U+ (X«)] + H.C. 

+ ^ [Ui (X^) [Uo (X«) X^U+ (X„)] + [Uo (X„) X^U+ (Xc,)] U+ (X„)] + H.C. - X« 

using now the fact that at the lowest order Uq (Xq,) Uq (X^)] = A^ (X^), one has : 
1 



= ^ [C^o (X«) [X«,t/o+ (X„)] + i/.C] 



1 

+ 2 



1 



aC/i (X,) (X,) + aAf (X„) t/+ (X,) + - [X,, U+ (X„)] + H.C. 



decomposing (Xa) into Hermitian and anti-Hermitian part we are thus led to : 

Aa-^ {Uo (X«) [X^,U+ (X„)] + H.C.) + [X^+aA^,ahr {U+ (X,))] 
+ f {A^hr {Ut (X,)) + hr {U^ (X„)) Af) 



(79) 



where ahr (Z) and /ir (Z) denote the anti-Hermitian and Hermitian part of an operator Z 
respectively. Now using Eq. (42), we are led to the following expressions : 



[X^+aA^,ahr{U+{X^))] = [B^,X^ + aAf-] 



and 



- {Afhr ([/+ (X„)) + hr {U+ (X„)) Af) + H.C. 

a {Afhr ([/+ (XJ) + hr (t/+ (X,)) Af) 
+ H.C. 



-A^ 



A^' A^' 
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with Ba given by the following relation 

B^ = - [.,£0 (X„)]-^ . (v. |^<"'Vi^„,£o (X«) + ^A''"Vp„,£o (X„) + H.C.^ 

{ [eo (X«) , - [^0 (X«) , + //.C.} j 
= - [., £0 (X«)]"' . (v^ {^<"Vp„,£o (X„) + ^<-Vp„,£o (X.) + //.C.j^ 

+ ^ |:p_^^"':p+^^"' - V-Aq-'v+Ao"' + h.c.^ (so) 

Note that given our previous notations, we can write: 

S„ = - [., £0 (X„)]-' . {(Px„£o (X«)) + H.C.} 

The two last expressions will allow to rewrite A^: To do so, we need to express carefully 
^ [Uo (X,) [X^,U+ (X,)] + H.C] . 

Actually, given our choice of symmetrization, at the second order this term is not equal 
to .4^ but also includes some second order corrections to symmetrize properly the product 
Uq (Xq) X [Xa,C/^ (Xq)] (see [10] for this fact). Moving powers of and on the left 
or on the right yields: 

^ {Uo (XJ [X^,U+ (X,)] + H.C.) 

^Af + ^A^-'Uo (X.) Vp„, ([X.,C/o+ (X.)]) + ^AtUo (X.) Vp„, {[X^,U^ (X.)]) 
= + I (A'^-t/o (X«) Vr^X + A^'^t^o (X„) Vp„,^^ + H.C.) 

-^A^"^ (Vp„,C/o (X.)) ([X„,C/o+ (X.)]) - (Vp„,C/o (X.)) {[X^,U^ (X.)]) + //.C. 

= Af + J (<"'Vp,,^^ + Vp,,^^ + H.C.) 

— - A^"^ A^"' A^ -\- - A^"^ A^"' A^ -\- H C 

= + ^ {At'^R^X + Vp„,^^ + Hc) + (^^Af [a^\A^] + H.c)j 
And as consequence, we can regroup the various terms of Eq. (79), to obtain ultimately : 

where we introduced the notations 

(r,+ ^X": Pa+f =A^" + J {^^' Vp„,^^ + A''"' Vp„,^^ + H.C.} 
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2. Hamiltonian diagonalization at the second order 



We can now turn to the expression of the diagonahzed Hamiltonian. Focusing first on 
the double integral term : 

h f-ai ^ 

'S'xr [Oq^5'x„^ [Oaa^o (Xaa)]] da2dai - - {Eq (X)) 



Jo 



./o 



a2 



V. 



+ 



(Px„,£o (X«J) Af"' + H.C^]] da^da, - ^ (sq (X)) 



^0 

~ {eo (X)) 



(Px„,£o (X«J) + H.C] ) ) + H.C.^ a^da 



due to the action of the shift operators, Xq^ and Xq,2 are shifted to X^ = X at that order. 

As a consequence, 

h />«! ^ 

'S'xft [Oai^x^^ [Oaj^o (XaJ]] (iaadai - - (^o (X)) 



Jo 



0.2 



h 



^0 



^ da^dar [V+ { (Px [V+ [(Px^o (X)) Af + H.C] ]) A^] + H.C] - - {so (X)) 

= y [n{ (^x [{T^^eo (X)) + ])A^}+ HC] - ^ {so (X)) 

We can now turn to the determination of the single integral term. For the sake of simplicity, 
we consider the case ((^ + (.)) H (X.a)) U^] — 0, which is in fact satisfied in the usual 
examples of interest.. 

/ 'S'xs [Oa£o (Xa)] da 

JO 

= ^ S^,V+ 1 1 {A^'VR,eo (X,) + Vr^So (X«) A^' + A^^Vp^e^ (X J + Vp.Eo (X„) 1 da 
+ I's^^n [^0 (X.) , A^] A:^ - \ [so (X.) , A^^] A^^ + HC^ da 

Introduce now the correction to the covariant derivatives at this order: T>^^ = (v^^jVp^^ 
with 
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Then, one can write : 



Jo 



'S'xa [Oa£o (Xa)] da= 



da 



which can be expanded as : 



'S'xfi [OaSo O^a)] da 



f 

Jo 

\ ^ [V+ [(Px„£o (X,)) Al'^ + H.C]] da 



V- I S^,V+i^[eoi:^a),jA^"'^x^,Af- + H.C.+ [B^,Ra + aAf"]\A^-+H.C.'jda 



S^,V+ ( [^0 (Xa) , Vx.,<" + H.C. + P„ + aA^"]] A^" + H.C^ 



da 



= -V^[{V^e^{X))Af + H.C] 



'-X'Vx,^^ + H.C.+ 



4 

<! ( so (X,) , Vx„,< + H.C. + 



I L (X) , ^A^'Vx^X + H.C. + 



B,l + A- 



Vx,£o {X) + H.C.] 
A^ + H.C. 
A^ + H.C. 



Px£o (X))^^ + //.C.] 



(81) 



with B = and 



h (\ 



AfVx.Ao'r 



4 V2 



(82) 



We also introduced the notation "Dx = ^"^r, "^p^ with 



Ultimately, one can gather the various terms to obtain the energy at the second order in h 
as a function of the canonical variables 



e (X) = eo (X) + \v. 



+y [n{ (^x [V+ [{V^e, (X)) Al + H.C] ]) A^] + HC] - ^ {e, (X)) 
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This formula is still compact and can be, for practical purposes, expanded in several ways. 
We first start by developing e (X) as a function of the zeroth order diagonalization data, 
that is Eq (X) and A^. Starting first by expanding A^, we are led to : 

e (X) = £o (X) + [V^SoAf + H.C] + y { (^x [ [V^SoA^ + H.C] ]) Af] + H.C] 



A^^VxX' + H.C.+ 



B,Ao"' 



e,,\A^^Vx^A^ + H.C. 



^0, 



A^ - 



Ar\ 



^0, 



1 



e,,-A^"'^x^,A^o+H.C. 



A^ + H.C. 



Af + H.C^ 



(83) 



Or, if we expand B fully, 

em 



eo (X) + ^V, 



{AfV^So (X)) + [eo (X) , X] X' - \ [^0 (X) , Al\ Af + B.C. 



f-2 

+-V+ { [AfV^V+ (a^V^So (X) + t [so (X) , + ff.C.)] + B.C.] 

+i-V+ { [V+ ((^^Vx£o (X)) + i [so (X) , ^^'] + B.C.) , } 
-ijV+ { [V+ ({AfV^So (X)) + i [e, (X) , A^^] A^^ + B.C.) , A^^] A^^ } 



A^ - 



Af + //.C. 



r r r r 



A^ + ii-.a 



■} 



(84) 



The last expansion is useful for practical purpose, since it yields directly the h expansion 
for the diagonalized Hamiltonian. However, it seems more elegant and relevant to rewrite 
e (X) as a function of the covariant variables x = (r, p) defined in the preceding section. 
Note first that the Berry phases Eq. 82, A^ = (^^ (X) , (X)) with A^ ^ Af + 



(j^'-VxiAf + [-B, X + A^] + B.C^ , satisfy, by construction : 



A"^' (X) 



-f 



S^^A^^da 



A^'{:^) = ljy^,A:^da 
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So that, recalling our previous definition of the variables x = X+^l-^^ (see Eq. 69), we can 
write for r. 



0<a<h 



0<ai<a 



+ H.C. 



daida 



R+hV+A'' (X) + - {{r+A^.Vn) V+A^ + {V+A^.Vp) V+A^ + H.C.) 



R+hV+Af + (Ia^'^xX + 



+ H.C. 



h 



- ( {V+A^.Vn) V+A^ + {V^AlV^) V+A^ + H.C.) 
= R+hA^ + 



(85) 



For p, we have in the same manner the following expansion 

p^P + J^^S^,V^Al^{R,P)da^+ J 5x,^ P+[^a-Vx J 7^+ K] 



0<a<?i 



0<ai<a 



+ H.C. 



daida 



P+hV+A'' (X) + - {{V+A^.Vn) P+< + {V+A^Vp) V+A^ + H.C.) 



+ H.C?j 



+ f (r+A^.Vn) V+A^ + (V+AlVp) P+< + H.c)j 



P+hA^ + ^AX 



As before, x = (r, p) satisfy of course an algebra Eq. (74). 
Now, come back to the compact form of e (X) : 



(86) 



e (X) = £0 (X) + ^V, 



(v^so (X)) + H.C. 



+f { (2^x [V^ [{V^eo (X)) A^ + H.C.] ]) A^} + H.C.] - ^ {so (X)) 
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expanded as : 



e (X) = £0 (X) + P+ I ^ (^^' Vx.eo (X) + V^.^o (X) ^^') | 

+^ [P+ { [eo (X) , A^^] - [eo (X) , A"^] A^^ ] + H.C] 

+-^7^+ { [AfV^V^ (^?Vx£o (X) + i [so (X) , ^^'] + //.C.)] + //.C.} 

+4^+ { ((-^^^x^o (X)) + I [so (X) , A'^'] X' + ^^-C^-) , A""'] X' ] 

To the second order, and given our choice of symmetrization, the contributions 

£o (X)+P+ 1^ (^^' Vx,£o (X) + Vx,£o (X) ^^') }+^^+ {^^Vxn (^^Vx£o (X)) + H.C.] 

can be recombined as : 

£o (^x+nn^""' + {{r+A^.Vn) V+A^ + {r+A^.Vp) V+A^ + i/.C.)^ = £o (x) 

On the other hand, one has also at the same order: 

^7^+ { [£o (X) , A^^] - [so (X) , ^^'] + //.a} 



+zy P+ { [AfV^r+ ( [^0 (X) , A^'] + i/.C.)] + H.C. } 



P+ ((^^Vx£o (X)) + i/.a) ,Af] A^^] 
-i^V+ { {{AfV^So (X)) + //.C.) , A^] A^'} 



7^4 



A^^ 



V+Af 



'-hV+ { [so (x) , i^'] i^' - [so (x) , i^'] i^' - [so (x) , [i^' , A'] ] } 



where we introduced the notations 



X' = \\^-\ (^+-4^ Vx)] (x) + H.C. 
Ultimately,the contributions to the Hamiltonian can be recombined to yield : 
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A' 



.a} 



|p+{[[eo(x),<' 



A^ 



Pi. 



So(x),»4q' Aq^Aq^ 

'so (x) , X^] , A'^'^ } - ^ (^0 (x)) (87) 

Note that in the all Hamiltonian, including the Berry connections A^^ and A^', we have 
replaced the operators (R, P) by x = (r, p) at each order of the expansion. 



VII. PHYSICAL APPLICATIONS 



A. The Dirac electron in an electric field 



To illustrate our general theory we consider the case of a Dirac electron in an external 
electric field. We will obtain the block diagonal Hamiltonian to the second order in h and 
will compare with the FW transformation obtained in [19]. Note that contrary to the FW 
which is not an expansion in h, the new method is valid for strong external fields (actually 
a FW transformation expanded into a power series in h was also recently proposed [20]). 

Let consider the following Dirac Hamiltonian (c — 1) 

Hi = a.P + ^m + eV (R) 

To compute the diagonalized Hamiltonian to the second order in h we first need the ze- 
roth order diagonalization transformation, which is the usual FW transformation for a free 
particle: 



E + ni + (3aP 



y/2E {E + m) 

where E — •\/P^~+~m^. In this case, we have (with S = 1 <t) 

-I3Q..VP + E{E + m) f3cx-iEP x E 



A^ = tUoVpU^' 



2£;2 {E + m) 



which leads to the first order projected Berry connections: 

P X E 



R 



2E {E + m) 



A^^V+A^^O 
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The zeroth order diagonalized energy is Eq (R, P) — E — -\/P^~+m^ 

Now, to complete the diagonahzation process, we need the matrix B Eq. (80) that enters 
in the definition of the covariant variables: 

Prom B and Eq.(82) we can then compute the Berry connections to the second order in 
h. Here expressions simplify greatly. Actually, one has : 

= A^ + ^[B,R+hA^] 

2E{E + m) 2 ^ 4^5 

t ^ -ci.VV + E{E + m)oi^ 
+^eVp ^^^^^ ^ — —^rV 



+he 



8E^ {E + m) 

E{E + m)Px {VrV X a) - (P X (P X a)) P.VrV 

4^4 {E + mf 



and 



A' 



we can also write the dynamical operators as: 



r= R+HV,A^ + -V, [B, A^] = + ^/^^" 

and 



4£;5 



Now using expression Eq. (87), we arrive at the following expression for the diagonal repre- 
sentation of the energy operator as a function of the covariant variables (r, p): 



A""' - L(x),i^'] i^'} 
= £o (x) + '-hV+ { [so (x) , ^^'] - [£o (x) , A'^q A^' } 
which, once developped, leads to: 

a = /3Vp2 + m2 + — Vr. + eV^ (r) 



(88) 



Here we have made the choice of fully symmetrizing in r and p, that is to weight equally all 
permutations in r and p in the series expansions of our expressions. One can now check that 
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developing the variables r as a function of the canonical variables R and P yields the same 
expression for the Hamiltonian as in [19] (apart from.a small sign mistake for the spin-orbit 
coupling in that reference) 

s = /JVP^ + eV (R) + ft^^|A_. 



8^5 



Note that this result can also be found directly by using the expression Eq. (84) for the 
diagonalized Hamiltonian as a function of the canonical variables. Eq. (88) being fully 
relativistic, to compare with the usual FW [15] approach we consider the non-relativistic 
limit and expand our result to second order in We readily obtain the well known diagonal 
representation of the positive energy (expressed in coordinates R and P) [15] : 

with cr the Pauli matrices. Note that it is the term order ft^ in Eq. (88) which, in the non 
relativistic limit, leads to the Darwin term. 



B. Massless Dirac particle in a static symmetric gravitational field 

We consider here the Hamiltonian of a massless neutrino propagating in an isotropic 
inhomogeneous curved space of metric ^''■' (R) = n"^ (R) 5*^ . (with the convention c—1) 

Ho^^ (a.PF(R) + F(R)a.P) (91) 

with F(R) = n-i (R). 

As for the electron in an electric field, we start by giving the zeroth order diagonalization 
matrix as well as the effective energy (R, P). Since at that order, (R, P) can be thought 
as commuting variables, one finds easily that the Hamiltonian diagonalization is performed 
through the following Foldy-Wouthuysen unitary matrix : 
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and 



So (P, R) = UoHo (P, R) U+ - (f(R)v^ + Vp^F(R)) 

We also need the Berry phases at the lowest order: 



(93) 



= I [t/oVpt/o+] 
= -i [C/VrC/+] = 



. -/3a.PP + £;2/3a - %EV x S 



The matrix B Eq. (80) needed to obtain the corrections to the energy at the second order 
is even simpler than for the electron in an electric field since here (see Eq. (82)): 



B^- [.,£0 (X.)]-^ .V. {^<"'V«„,5o (X) + ^V«„,£o (XX^'I = 

Therefore Berry connections and covariant variables to the second order reduce to 
and = as well as r= B^^KP^J^ = K+h^^, P = P- 

Now, using expression Eq. (87), we ultimately obtain the energy operator: 

£ (P, r) = 1/? (f (r) Vp~^ + Vp~^F (r)) - ^ (^o (X)) 



•^0 



(^F(r) Vp2 + A/p2F(r)) - ^/i \/p2, F (R) 
(f (r) Vp^+ Vp^F (r)) - ^^PVF (r) 



(94) 



4 a/P2 

From the last formula, we can deduce the equations of motion to the second order approxima- 
tion. Indeed, for a particle of positive energy only e (P, r) = | (f (r) VP2 + Vp^F (r)) - 



1 K-" 



^^_^PVF(r) with now r = R+i^^ a 2 x 2 matrix (the spin matrix is S = her/ 2), the 
usual relations r = — | [r, £ (P, r)] and P = — [P, £ (P, r)], lead to the equations of motion 



r = Vp£ + /iPxe' 

P = Vr£ 



(95) 



Note that the anomalous velocity term hPxQ^'^ in Eq. (95) has to be understood as a 
symmetrized expression of P and r. Here, [ri,rj] = ifi^9[J = ifi?eijkQ]I = —ifi^SijkX^, and 
A = cr.P/P the helicity. In addition helicity is conserved d\/dt = 0. Eqs (95) imply the 
spin Hall effect of light (for the photon one just has to replace the Pauh matrices cr by the 
spin one matrices) as a consequence of the anomalous velocity term KPxQ^^ (see also [26] 
and [27] for the experimental confirmation of this effect) 
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Prom the Eq. (95) we deduce the following expression for the velocity components 

. _ 1 /_^P^ P^_^\ Xh_ / fe glnn c dlnn c A 
~ 2 \n{r) P n(r) ) ^ 2P^^'''' \ dx' n{r) ^ Ox' n(r) ) 

from which we compute the modulus of the velocity 

^ = ;^(l + ^'^((Vlnn)^--l(P.Vlnnf)) ^ + Oih') (97) 

This expression shows a very small correction to the usual expression v = c/n{r) of order 
due to the interaction between the polarization state and external inhomogeneities. 

VIII. CONCLUSION 

In this paper, we presented a new diagonalization method for a generic matrix valued 
Hamiltonian which leads to a diagonal representation where the operator energy takes an 
elegant and compact form. This approach requires the introduction of some new mathe- 
matical objects like non-commuting operators evolving with the Planck constant promoted 
as a running variable and thus reveals a mathematical structure reminiscent of the stochas- 
tic calculus. It also shows once more the very important role played by Berry phases in 
these systems as the energy operator is written in terms of covariant dynamical operators 
containing Berry connections and satisfying a non-commutative algebra. 

It was also found that the diagonal representation of the energy is solution of a differential 
equation in h presented previously in [8] and which could only be solved recursively in a 
series expansion in h. Actually the formal exact solution presented here can also be written 
explicitly as a scries expansion in H, but it appears that the derivation of the coefficients 
of this expansion is now much more easier. Indeed we could give the expression of the 
energy and the dynamical variables to the second order in h for a generic matrix valued 
Hamiltonian. We then apphed this method, first to the simple case of a Dirac electron in 
an external electric field which allowed to recover the usual Pauli Hamiltonian in the non 
relativistic limit including the Darwin term of order h"^, and second to the neutrino in a 
gravitational field. This is obviously a good check for the validity of the proposed method. 
We leave for subsequent work its application to more complicated systems in condensed 
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matter or relativistic particle physics such as Bloch and Dirac electrons in interaction. 
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